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Abstract 

We prove estimates for charge transfer Hamiltonians, including matrix and inhomoge- 
neous generalizations; such equations appear naturally in the study of multi-soliton systems. 

1 Introduction 

This paper is devoted to the study of dispersive properties of the model corresponding to the 
time-dependent charge transfer Hamiltonian 



m 

H{t) = --A + J2v,{x-Vjt) 



with rapidly decaying smooth potentials Vj{x) and a set of mutually non-parallel constant velocities 
Vj. Our main focus is on the decay estimates for the solutions of the time-dependent problem 

-dtiJ + H{t)'il^ = 

associated with a charge transfer Hamiltonian H{t). The well-known estimates for the free 
Schrodinger equation {H{t) = — ^A) on M" are 

l|e**^/||LP<Cp|tr"^^"^^||/|li.', P>2, - + ^ = 1. 

p p 

They imply the Strichartz estimates 

\\e''-fh^^Lr<C,\\fh2, 2<r,g<oo, - + - = -, n>3 [GV,KT] 
t X r q I 

Such estimates play a fundamental role, among other things in the theory of nonlinear dispersive 
equations. The extension of such theories to inhomogeneous problems (either due to curvature, 
local potentials or coherent structure such as solitons, vortices etc. ...) then motivated the efforts 
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to establish the decay estimates for the general time independent Schrodinger operators of the 
type H = — + V{x). In this case the estimates should take the form 

(1.1) lle-'^^PcWV'ollLP <Cp|tr"(^-i)||V^o|lLp' P>2, - + ^ = 1, 

p p' 

(1.2) ||e-^*^Pe(^)/||r«L. < QII/IIl^ for 2 < < oo, - + - = 5, n > 3 

I ^ r q z 

where Pc{H) is the projection onto the continuous part of the spectrum of the self-adjoint operator 
H. Its purpose is to prevent the emergence of any bound state, i.e., an eigenfunctions of H. 
Under the evolution e~^^^ such bound states are merely multiplied by oscillating factors and thus 
do not disperse. 

The first approach to the proof of estimates was developed by Journe, Soffer, Sogge |JSSj . 
They used a time dependent method which combined spectral and scattering theory with harmonic 
analysis. Their method involved splitting solutions into high and low energy parts and using 
Kato's smoothing and the local energy decay on the corresponding pieces. Later, a stationary 
method was used by Yajima |Yalj who proved that the wave operators are L^-bounded under 
weaker assumptions on the potential than in |JSSj . His theorem implies the dispersive bounds. See 
also Weder |Wej for results in one dimension n = 1 and Yajima |Ya3j for n = 2. 

Most recently, Rodnianski and Schlag |RSj have addressed the issue of determining the optimal 
class of potentials V{x) for which the estimates (ll.lj) and (|1.2j) hold true. In particular, they 
obtained L°° decay estimates for the scaling invariant "L^'-like" classes of small potentials in 

dimension n = 3. In |NSj . Nier and Soffer were able to establish the decay and Strichartz estimates 
for finite rank perturbations of Hq = — y • The Strichartz estimates for the inverse-square potential 
were obtained by Burq, Planchon, Stalker, and Tahvildar-Zadeh in |BPSTj . 

The situation in the case of time-dependent potentials V{t^ x) becomes more complicated. On 
a superficial level, this can be attributed to the disappearance of a connection between the decay 
estimates and spectral theory of a corresponding Hamiltonian. In particular, solutions that remain 
trapped in a compact region, and thus do not disperse, might still exist. However, they can no 
longer be easily characterized as corresponding to the bound states of some Hamiltonian as in the 
case of time independent Hamiltonians. 

Let U (t) be the solution operator corresponding to the time-dependent Schrodinger equation 

(1.3) -dt^-^iP + V{t,x)iP = 0, 

I 2 

V'|t=o = V'O, 

i.e., ip{t, •) = U{t)ipo. Then the desired decay estimates take the form 

(1.4) r(t)^o||Lp <Cp|tr"(^-i)||Vo|li/ P>2, - + ^ = 1 

p p 

holding for an appropriate set of initial data. The first results in this direction were obtained 
in |RSj . where the L°° decay estimates for all initial data were established in dimension 

n = 3 for a large class of small time-dependent potentials. The smallness assumption rules out the 
emergence of trapped solutions. 
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In this paper we prove decay estimates for solutions of a time-dependent Schrodinger equation 
corresponding to a charge transfer Hamiltonian 

H{t) = -^A + ^ Vj{x - Vjt), Vj j^Vi, i / j, 
j 

see Theorem 12.71 below. Our main motivation comes from the problem of asymptotic stability of 
noninteracting multi-soliton states. This question refers to solving a NLS 

(1.5) idt^P + + = 

in M", n > 3 with initial data tpo = Ylj=i "^ji^^ ') where Wj are special standing wave solutions 
called solitons and Rq is a small perturbation. In our forthcoming paper |E,SSj we show that if the 
solitons are sufficiently separated at time t = 0, and if Rq is sufficiently small in a suitable norm, 
then the solution ip evolves like a sum of solitons with time-dependent parameters approaching a 
limit, plus a radiation term that goes to zero in L°°(R^). This argument requires linearizing the 
NLS around the bulk-term, which is given by the sum of the solitons. This leads to the problem of 
establishing dispersive estimates, which typically means L°° decay estimates, for the linear 

problem. It is well-known that the linearized equation needs to be written as a system in 5iJj and Sip 
where Sip is the variation of the bulk-term. A first step in the understanding of this problem is to 
consider the scalar model that is closest to the system at hand. It is easy to see that these scalar 
equations are Schrodinger equations of the charge-transfer type 

(1.6) -di; --,!; + Y,Vj{--tv,)i; = 0, 

' ^ j=i 

with real potentials Vj and distinct vectors Vj G M". 

This problem has been extensively studied in the literature in connection with the question of 
asymptotic completeness. Asymptotic completeness for the equation (|1.6j) had been first established 
by Yajima |Ya2j . Later, Graf |Grj has obtained a new proof containing a crucial argument showing 
that the energy^ associated with a solution of (|1.6|1 remains bounded in time. For similar results 
see also the work of Wiiller |Wuj . and for more general charge transfer models with PDO's, the 
work of Zielinski [Zij . Their results imply that for a dense set of initial data, solutions converge to 
the states represented by radiation and the sum of the bound states, associated with each of the 
Hamiltonians Hj = + Vj{x), traveling with velocity Vj. 

In this paper we obtain decay estimates for (|1.6|) . which in particular lead to a very short 
and self-contained proof of the asymptotic completeness theorem. The idea is to decompose the 
solution into k+1 channels that move (and expand) along with the potentials according to a straight 
motion. On the j^^ channel one has well-known dispersive estimates by Journe, Soffer, Sogge |JSSj 
and Yajima |Yalj for the equation with the single potential ^ — vjt), which dominates on that 

^Energy here simply means H^fi for a solution While the conservation of the norm is a trivial 

consequence of the potentials Vj being real, the energy bound is by no means a simple matter since the problem 
contains a time-dependent potential. 

^The results of |JSS| and |Yall apply to problems with time-independent potentials. However, the problem with 
a single potential Vj{x — Vjt) can be reduced to a problem with the potential Vj{x) by means of a Galilei transform. 
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channel. One then proceeds to show that the interaction between the channels is weak enough not 
to destroy this property. This weak interaction is accomplished by distinguishing low from high 
momenta. In the former case one uses the slow propagation rate, in the latter a variant of Kato's 
smoothing bound, to derive the desired bounds. As always, we cannot hope to prove dispersive 
bounds for all initial conditions. In the scalar equation with a single stationary potential this simply 
means being perpendicular to the bound states. In the charge transfer case one expects this to be 
replaced with perpendicularity to the (traveling) bound states, or at least that the projection of 
the solution onto the traveling bound states should go to zero. This is precisely the "asymptotic 
orthogonality condition" which we use below. 

For this type of initial data we establish the decay estimates 

(1-7) \m)\\L^<C{n)\t\-'^\\M\L^nL^, 

which, by interpolation with a trivial bound, also imply the the full range of the estimates 
(1.8) \\m\\Lv<C{n,p)\t\~'^^^-'v'^Uo\\Lv'nL^^ 2 < p < oo 

The presence of the LP' norm on the right-hand side of (|1.7|) and (|1.8|) distinguishes them from 
the standard — > estimates. This, however, is no more than a technical issue. In fact, our 
main result is stated as an n ^ + L°° estimate, which, for example, is sufficient for a 
proof of asymptotic completeness. We then establish a procedure converting this estimate into an 
n — > L°° bounds, using an argument along the lines of |JSSj . 

To motivate the results of the second part of this paper, we recall that stability and asymptotic 
stability of solitons and other coherent structures require analyzing the spectral properties, scat- 
tering, and local decay of the linearized matrix operator. A comprehensive description of results 
in this direction can be found in the book by C. Sulem and P. L. Sulem |SuSuj . In particular, the 
spectral properties of the linearized operators for NLS solitons were developed in Weinstein |Weinj , 
Shatah, Strauss |ShStj . and Grillakis, Shatah, Strauss [GrShStj . The linearized theory of other 
structures was also considered, e.g., vortices in Weinstein, Xin [WeinXij . and Hartree-type solitons 
in Frohlich, Tsai, Yau |FTYj . Many of these results are examples to which our LP theory applies. 

On a more technical level, the second part of this paper extends our approach to the full system 
that one obtains by linearizing H1.5() around a sum of solitons. Several complications arise, even for 
a system with a single stationary matrix potential: Firstly, these systems are no longer self-adjoint 
and there is the well-known issue of existence of unstable modes, see Weinstein |Weinj . Secondly, 
and somewhat related to the aforementioned instabilities, the L°° decay estimates are more 

involved for systems ( j.TSSj and (Yal) do not apply). 

As far as the latter is concerned, Cuccagna [Cuj recently established the necessary bounds by 
transferring Yajima's L^-wave operator approach to the system case. While the wave operators 
are of course of independent interest, we feel that it is desirable to have a more direct approach to 
the dispersive estimates for systems. It turns out that one can easily adapt a simple but powerful 
method of Ranch |Q to systems, which is what we do in Section[7|below, see Theorem 17.41 Ranch's 
method, which is related to earlier work of Dolph, McLeod, Thoe |DMcLT] and Vainberg |Vaij . 
hinges on analytic continuation of the resolvent across the spectrum as a bounded operator in an 
exponentially weighted L^-space. This can only be done for exponentially decreasing potentials, 
but the potentials arising from linearizing around solitons are of this type (since solitons decay 
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exponentially). It is important to realize, however, that Ranch's method only leads to local L?- 
decay, which is inadequate for the nonlinear applications. To obtain L°°-bounds, we then use an 
observation of Ginibre |Ginj that allows one to pass from local L^-decay to L°°-decay by means of 
a double application of Duhamel's formula. 

After deriving these bounds for systems with a single stationary matrix potentials, we then 
introduce charge transfer models in the system case, see Definition 18.11 These are of course anal- 
ogous to the scalar models ()1.6() but differ in that one not only needs to translate the potentials 
according to linear motion, but one also needs to modulate the matrix potentials in a suitable 
manner. In particular, these matrix potentials have complex off-diagonal entries. We then proceed 
to apply the aforementioned method of "channel-decomposition" to these matrix charge transfer 
equations, which yields bounds of the form (|1.7() . Finally, we extend the decay estimates to the 
inhomogeneous equation with charge transfer Hamiltonians. We give a general L^-decay estimate 
in terms of a properly chosen norm of the inhomogeneous term, see Theorem 16.11 

2 Charge transfer model 

Definition 2.1. By a charge transfer model we mean a Schrddinger equation 

(2.1) -dt^P - l-A^j + K(x - v^t)^|J = 

^ K=l 

where are distinct vectors in M", n > 3, and the real potentials are such that for every 
1 < K < m, 

1. has compact support (or fast decay), Vk,Wk G 

2. is neither a zero eigenvalue nor a zero resonance of the operators 

H^ = -^A + V^{x). 

The regularity assumption on the potential can be relaxed considerably without too much effort, 
but we use the C^-property for convenience. The assumption n > 3 will be made throughout. This 
definition is standard, see |Grj . |Ya2j . The Schrodinger group e"**^" is known to satisfy the decay 
estimates (see Journe, Soffer, Sogge |JSSj and Yajima |Yalj ) 

(2.2) \\e-''''-P,{HM\L^ < ir"/'ll^o||Li. 

Here Pc{Hi^) is the spectral projection onto the continuous spectrum of H^. The estimate (|2.2|) 
holds under the "no eigenvalue/resonance" assumption 2. above for potentials obeying conditions 
roughly equivalent to the assumption 1. 

Remark 2.2. It follows from the results of Yajima |Yalj that in dimension n = 3 the estimate 1)2. 2(1 
holds for any bounded potential with sufficiently fast decay at infinity. In higher dimensions one 
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needs to add assumptions on the derivatives of the potential V,^: 

(2.3) |V°y,(x)| < C„(x)-^ v|a|<^^-3, S > max{n + 2,^ - 2}, 

n-i 2(n-2) 
n — 1 



We shall assume that potentials Vk{x) always obey conditions guaranteeing estimate (|2.2j) . e.g. 
conditions (fTl)) . 

Definition 2.3. Xei + L°° := {/ : 3h; g : f = h + g, \\h\\i^2 + ||(7||loo < oo} with norm 

■= inf/=h+g(||/i||L2 + IIsIIl-) 



We shall use a weaker version of the decay estimate which follows from ()2.2p combined with the 
unitarity of e~**^'*, namely 



(2.5) ||e-^*^«Pe(F,)Vo||L-+L2 < {ty^'m 



LinL2 



where {t) = (1 + t^)2 . In the following we shall assume that the number of potentials is m = 2 and 
that the velocities are vi = 0,V2 = (1,0, ... 0) = ei. This can be done without loss of generality. 
An indispensable tool in the study of the charge transfer model are the Galilei transforms 



-i^-^t —ix-v i{y+tv)-p 



(2.6) Q^^yit) = e-'—'e 
cf. |Grj . where p = — iV. These are the quantum analogues of the classical Galilei transforms 

x>-^x — tv — y, p^p — v 

in the following sense: if / is a Schwartz function, say, such that / and / are supported around 0, 
then Qij^y{t)f is supported around tv + y, and Q^^y{t)f is supported around v. Under 0^,y(t) the 
Schrodinger equation transforms as follows: 

Qv,y{ty^^ = e**^S^,;;(0) 

and moreover, with H = — + V , 

(2.7) m := 0^;,j,(t)-'e-**^0,^,j,(O)(/)o, , S.-y(i)~' = e~'yU-v,-y{t) 
solves 

(2.8) -dtil^ - ^ + V{- -tv-y)i) = Q 



These properties are not only easy to check, but are of course to be expected in view of the 
classical interpretation. We will make frequent use of these transformation laws without further 
notice. Another property that we will use often without further mention is that the transformations 
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Qv,y{t) are isometries on all spaces. Finally, since in our case always y = 0, we set 0i;(t) := 0i;,o(O- 
By (|2.7|) . 0ei(*)~^ = 9~e{{t) in that case. 
We now return to the problem 

(2.9) -dt'ilj-lA'4; + Vi'ifj + V2{- -te{)i; = 

with Vi,V2 compactly supported potentials. Let ui,...,Um and wi,...,W£ be the normalized 
bound states of Hi and H2 corresponding to the negative eigenvalues Ai, . . . , Am and fii, . . . , m, 
respectively (recall that we are assuming that is not an eigenvalue). We denote by P}j{Hi) and 
Ph{H2) the corresponding projections onto the bound states of Hi and H2., respectively, and let 
Pc{Hk) = Id — Ph{Hfi), K = 1,2. The projections P;,(-f^i,2) have the form 

m t 

Pb{Hi) = '^{■,Ui)ui, Pb{H2) = y^j-, Wj)wj. 
i=i j=i 

In order to state our main theorem, we need to impose an orthogonality condition in the context 
of the charge transfer Hamiltonian (|2.9)) . 

Definition 2.4. Let U{t)ipQ = il){t,x) he the solutions of (|2.9)) . We say that ipQ (or also ^/'(t,-)j is 
asymptotically orthogonal to the bounds states of Hi and H2 if 

(2.10) \\PbiHi)U{t)^Po\\L2 + \\PbiH2,t)Uit)ML2 ^ as t ^ +00. 
Here 

(2.11) Pb{H2,t) := Q_^^^{t)Pb{H2) Qeiit) 
for all times t. 

Remark 2.5. Clearly, Pb{H2,t) is again an orthogonal projection for every t. It gives the projection 
onto the bound states of H2 that have been translated to the position of the potential V2(- — tel). 
Equivalently, one can think of it as translating the solution of 1)2. 9() from that position to the origin, 
projection onto the bound states of H2, and then translating back. The explicit form of Pb{H2,t) 
is given by 

^ r 

(Pfe(/72,t)/)(x) = Ve*"i«;,(x-tel) / fiy)e-'y'wjiy-te-'i)dy. 
j=i 

We will make little use of this formula, though. 

Remark 2.6. It is clear that all ipQ that satisfy (|2.1U|) form a closed subspace. This subspace coincides 
with the space of scattering states for the charge transfer problem. The latter is well-defined by 
Graf's asymptotic completeness result jHll- This is discussed in more detail in Section^ 

We now formulate our main result. 
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Theorem 2.7. Consider the charge transfer model as in Definition \2. 1\ with two potentials, cf. H2.9() . 
Let U{t) denote the propagator of the equation ()2.9|) . Then for any initial data ipo S L^nL^, which 
is asymptotically orthogonal to the bound states of Hi and H2 in the sense of Definition \2.4\ one 
has the decay estimates 

(2.12) \\U{t)ML^^L^ < W/'ll^ollLinL^- 

An analogous statement holds for any number of potentials, i.e., with arbitrary m in ()2.1|) . 

We prove this theorem for the case of two potentials, but this is for simplicity only. The same 
argument also apphes to the general case of m > 2 (one then needs to split into m + 1 channels, 
see below). A more substantial comment has to do with the assumption of compact support of the 
potentials. Inspection of the argument in the following section reveals that it equally well applies to 
exponentially decaying potentials, say. But also power decay is allowed, provided it is sufficiently 
fast. We shall prove H2.12|) by means of a bootstrap argument. More precisely, we prove that the 
bootstrap assumption 

(2.13) \\U{t)M\L^+L^ < Co(t)-"/'||V'o||LinL2 for ah < i < T 
implies that 

(2.14) mmoh^+L^ < Y(*)""^'llV'ollLinL2 for ah 0<t<T. 

Here Co is some sufficiently large positive constant which is assumed to be much bigger than any 
universal constant appearing in our calculations (e.g. ||Vi,2||LinL°° > implicit constants in Propo- 
sition |S3 below, and decay bounds (|2.5|) 1. The logic here is that for arbitrary but fixed T, the 
assumption 1)2. can be made to hold for some Co, depending on T. Iterating the implication 
(|2.1H|) =^ (|2.14|) then yields a constant that does not depend on T, thus proving the theorem. We 
shall continue to use the notation < to denote bounds involving multiplicative constants indepen- 
dent of the constant Co. 

Remark 2.8. Since 

\\U{t)i^o\\L2+L°- < \\U{t)'4)o\\L2 = \\iPo\\l^ < Wi^oh^nL^ 
it suffices to prove (jTTll) for t > to := {^ f^"' . 

Remark 2.9. For the nonlinear applications it might be also useful to have estimates as in Theo- 
rem EIZI for perturbed charge transfer Hamiltonians. It is easy to see that the method of proof is 
stable under such small perturbations, see Section [3.71 below. 

3 Proof of the decay estimates 
3.1 Bound states 

Our first result concerns the rate of convergence of the projections onto the bound states of solutions 
which are asymptotically orthogonal to the bound states. 
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Proposition 3.1. Let ij^it^x) = {U{t)'ipQ)(x) be a solution of H2.9|) which is asymptotically orthog- 
onal to the bound states of Hi and H2 in the sense of Definition\2.^ Then 



\\Pk{Hi)U{t)M\L^ + \\Pb{H2,t)U{t)M\L^ < e-"1^o||L2 

for some a > 0. 

Proof. By symmetry it suffices to prove the bound on the first part. More precisely, let U{t) := 
Sei(i)f/(t), and (t){t) = U{t)%l)o. Then (t){t) solves 

(3.1) -dt<t) - J A,/. + Vi{- + te{)(k + V2(t> = 0, 

0|t=o(a;) = (Sei(0)V'o)(a;), 

and 

||Pfe(i?2,t)t/(t)V'||2 = ||P6(//2)0el(t)f^(i)V'o||2 = ||n(^2) V'O || 2 • 

This clearly allows one to reduce the treatment of to that of fl\ . Decompose 

m 

(3.2) \J(t:)i)Q = ^a,{t)ui^^i{t,x) 

1=1 

relative to Hi so that 'ipi{t,-) lies in the continuous subspace of Hi, i.e., Pc{Hi)iJji = ipi and 
Pb{Hi)ipi = 0. By assumption, 



|aj(t)|^ ^ as t — > DO. 



1=1 



Substituting (glS) into yields 



-dtil^i - JaV'i + Vi-^i + V2{- - te"I)V'i+ 
I 2 



(3.3) + 



^1 , , 1 



-dj{t)uj — -Aujaj{t) + Vi Ujaj{t) + V2(- — tei)ujaj{t) 



0. 



Since Pc{Hi)ipi = ipi one has 

(_1a + Vi)yJi = Hi^i = P,{Hi)Hi^i, dtiPi = Pc{Hi)dtil>i. 

In particular 

PbiHi) - \Mi + ^iV'i) = 0. 

Thus taking an inner product of the equation with and using the fact that {uK,Uj) = 6j 
as well as the identity 

--Auj + ViUj = XjUj 

we obtain the ODE 

^ m 

-a,,{t) + Xna^,{t) + {V2{- - tei)7pi,u^,) + 'Yaj(t){V2{- - tei)uj,u^,) = 
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for each with the condition that 

a^it) — > as t ^ +00. 

Recall that is an eigenfunction of Hi = —^A + Vi with eigenvalue < 0. It is well-known (see 
e.g. Agmon |Ag| ) that such eigenfunctions are exponentially localized, i.e., 

(3.4) / e^^l^l \uK{x)f dx < C = C(Vi,n) < 00 for some positive a. 
Therefore, 

(3.5) ||'^"2(--te"I)'u^||2<e-"* for all t>0 

which follows from the assumption that V2 has compact support. The implicit constant in 1)3. 5() 
depends on the size of the support of V2 and ||1^2||l°°- Therefore, 

Ut) := {V2{- -te-{)xPi,u^) 

satisfies 

(3.6) \Ut)\ < e-^'WMi^ < e-^'\\P,{Hi)U{t)ML^ < e-'IV'olb- 
In view of ()3.1|) . solves the equation 

^ m 

(3.7) -a,(i) + X^a^it) + ^ aj(t)Cj«(t) + f^{t) = 

' i=i 

0^(00) = 0, 



where Cj^it) = C^j{t) = {V2i—tei) ). By maXj fj \CjK{t)\ < e~°*. Solving (jnHj) explicitly 
we obtain 

a{t) = ie-' -fo ^(^) '^^ / -fo ^(^) <^^f{s) ds, 



where Bj^{t) = Xj6j^ + Cj«.(t). Hence, by unitarity of e* i'o -^(^^ "'^ and (ESJ, we conclude that 

|/(s)|d.<a-ie-"*||V'o||L2. 
This proves the proposition. □ 
3.2 The three channels 

Fix some small S > 0. We introduce a partition of unity associated with the sets 
^5^(0) = {x : \x\ < 5t}, Bstitei) = {x : \x - tei\ < 6t} 

and 

R''\{Bst{0)UBst{te-{)). 
10 



Namely, let xiit^^) be a cut-off function s.t. 

Xi{t,x) = l, xeBstiO) and xi(.t,x) = 0, x e R"" \ B2St{0) ■ 
Here t > to ^-nd 5 > is a fixed small constant. Define 

X2{t,x) = xi{t,x- tei), X3{t,x) = 1 - xi{t,x) - X2{t,x). 

Observe that since to = (^)^''" is large, see Remark [2.8l the support oixi{t, •) contains the support 
of Vi for all t > to and the support of X2{t, •) contains the support of V2(- — tel). 
There is the following natural decomposition of the solution U{t)ipQ: 

= Xiit, ■)Pb(,Hi)U{t)i,o + Xiit, ■)Pc{Hi)U{t)i;o + X2{t, •)n(^2, t)U{t)i;o 
+X2it, ■)Pc{H2,t)U{t)iPo + X3{t, ■)U{t)i^o- 

The terms xiPbiHi)U{t)ipQ and X2Pb{H2,t)U{t)'ipQ can be estimated immediately using Proposi- 
tion 13.11 Since U{t)ijjo is asymptotically orthogonal to the bound states of both Hi and H2 by 
assumption, one concludes that 

||xi(tr)n(^^i)^7(t)^o||L2+L- + IIX2(tr)n(^2,t)i7(t)V'o||L2+ioo < e-^^llV'ollLinL^. 

To estimate the remaining three terms 

Xi{t,-)Pc{Hi)U{t)tPo, X2{t,-)Pc{H2,t)Umo, X3{t,-)Umo 

we use Duhamel's formula. However, on each of the three "channels" we compare U{t)ipQ to a 
different group. We claim that on the support of Xi{t-, ) the solution U{t)ipQ is best approximated 
by e~**^iVo) on the support of X2{t, •) by e'^'^^^ifjo, and on Xsit^ ') by e^*~'ipo. More precisely, with 
Xi = Xi{t,x), 

(3.8) xiPc{Hi)U{t)i;o = Xie-'^^^PcC^OV-o - iXi f e-^(*-^)^iPe(^i)^2(- - se-{)U{s)ijQ ds. 

Jo 

To analyze the second channel, define U{t) := (t)C/(t) g_e|(0). Then cj}{t) := U{t)(j)o solves 
equation (|3.1() . where (j)Q := 5el(0)V'o- The Duhamel formula relative to e~^^^'^ can now be written 
in two different ways, namely with respect to a stationary frame (this involves U{t)) or with respect 
to the moving frame (this involves U{t)). Indeed, 

(3.9) xiPc(^2)f/(t)0o = Xie-**^^Pc(^2)</'o - iXi f e-"^'-'^""^ P,{H2)Vi{- + sel)C/(s)(/>o ds 

Jo 

X2Pc{H2,t)Uit)i;o = X2 0-ei(t)e-**^^Pc(i^2)Sei(O)Vo 

(3.10) -iX2 9-el{t) [ e-'^'~'^''^Pc{H2)Vi{- + se-l)Q,-i{s)U{s)i^ods. 

Jo 

The stationary formulation (|3.9() is basically the same as 1)3. 8() . whereas (|3.10() follows from 1)3. 9() 
by means of (|2.11j) from Definition 12.41 It is evident from (|3.8|) and (|3.9j) that it suffices to prove 
the requisite amount of decay for xiPc{Hi)U{t)'ipo. However, we will need p.lUj) in the treatment 
of the first channel, i.e., (|3.8|1 . Finally, for the third (the "free") channel, one has 

(3.11) X3Umo = X3e''^^Po-iX3 f e'^'-''^^ {Vi + V2{- - se-{))U{s)iPo ds. 

Jo 
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3.3 Analysis of xi/^c(^i)^^(t)^o 
We start with the following simple lemma. 

Lemma 3.2. Let W £ LP n L^ for some p G [1, 2]. Then 

\\Wf\\Lr'<\\W 



Ln particular, the dual space of L + L°° is L H L . 

Fix sufficiently large constants A, B > 0, B >> A, which are independent of the constant Co 
and write the following expansion 



XiPc{Hi)U{t)i;o = Xie-**^^Pc(^i)V'o-iXi 
(3.12) 

According to (|2.5j) 



t-A 



e-*(*-^)^ip,(i7i)y2(- - se-{)Uis)iJods 



iXi / e-*(*-^)^iPe(^i)V2(- - sel)C/(s)Vo ds. 



(3.13) 



Using 1)2. 5() . Lemma 13.21 and the bootstrap assumption (|2.13p 



t-A 



e-*(*-^)^iPc(^i)V2(- - sel)C/(s)Vods 



< 



< 



t-A 



- 



t-A 



- s)'^/^ 



\\V2{- - sel)C/(s)V'o||LinL2 
l|V2||LinL2 + ll^2||L2nL-' \\U{s)4)o\ 



ds 



t-A 



{t - s)"/2 (s)W2 



ll^2||LinL°°ll'0o||LinL2 



+ 



1 



(3.14) 



ds\\V2\\LinL^\\Uis)i;oh 

+ A{ty 



/o 

<Co^-("/2-i)(t)-"/2||y2||LinL^llV'o||LinL2 



^0112- 



Thus, since n > 3, there exists a choice of a sufficiently large constant A that does not depend on 
Co such that 



(3.15) 



t-A 



e-^(*-^)^iPe(ii'i)V2(--sel)?7(s)V'ods , < 10-2Co(t)-"/2||^o||^,^^, 



L2+L° 



provided Co was chosen sufficiently large. Recall that t > to = (Co/2)2/". We can therefore 
assume that t >> A. The third term on the right hand-side of 1)3. 12() is treated as follows. With 
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f e-'^'-'^"^Pc{Hi)V2{- - se{)U{s)iPo ds 
Jt~A 



Xi 



Jt-A 

(3.16) +xi r e-^(*-^)^ip,(i/i)y2(- - se-{)P,{H2, s)U{s)4^o ds =: + J,. 

Jt-A 

Since U{s)ipo is asymptotically orthogonal to the bound states of both Hi and H2 by assumption, 
Proposition 13.11 imphes that 

\\PbiH2,s)Uis)^Po\\L^ < e-'^^WMi^- 

Therefore, using unitarity of e^*(*^*)^i and boundedness of the spectral projection Pc{Hi) we 
obtain 

WJbh^+Loo < \\M\L^<f \\V2{--se-{)Pi{H2,s)U{s)ML^ds 

Jt~A 

< \\V2\\l^ f \\Pi{H2,s)U{s)i>Q\\L2ds 

(3.17) < \\V2\\l^ f e-'''dsUo\\L^<\\V2\\L^Ae-'''/^\\M\L^nL^, 

Jt-A 

assuming as we may that t — A > t/2. This imphes the bound 

(3.18) llJfellL^+L- < (i)""/'llV'o||LinL- 

To deal with the term Jc in H3.16() we expand X2{s, •)Pc{H2, s)U{s)^po via the Duhamel formula H3.1U() . 
Then 

Jc = XI t e-^(*-^)^ip,(//i)y2(--sel)0,,l(s)e-^^^2P,(if2)0ei(O)V'ods 
Jt-A 

-iXi f re-^(*'^)^^Pc(^i)^2(--sel)0_ei(5)e-*(^-^)^^P,(i/2) 

Jt-A Jo 

Vi{- + rel) 0el {t)U (r)V'o drds 

(3.19) =: + Jl 

Using (|2.5() for k = 2 and Lemma 13.21 we estimate as fohows: 

\\Jc\\l^+l^ < \\Jc\\L^<f V2{--se-{)Q.,-',{s)e-''''^P,{H2)Qe-Mi'ods 

Jt-A 

< imi^nLoo f ||0-el(5)e-'^^^Pc(i?2)0el(O)Vo||L2+ioods 
Jt-A 

(3.20) < \\V2\\L^nL^ j'_^j^\mL^nL^ ds < \\V2\\L2nLo^A{t)-^/^i^o\ 



LinL2. 
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To pass to (|3.20|) one uses that 

\\Pc{H2)f\\L^nL^ < ll/llLinL^ + mmfh^nL^ < Wfh^nL^. 
Since ^ is a fixed constant independent of tlie constant Co we conclude tfiat 

(3.21) ll4llLHL-<(t)-"/'llV'o||LinL- 
We decompose as follows: 

Jt-A Jo 

Vi{- + rel) gei(r)?7(r)V'o drds 
+Xi{tr) f f e-^(*-^)^ip,(^i)F(|p-1<M)y2(--5el) 

Jt-A Js-B 

0_el(5)e"*^'~"^'''i^c(i^2)l^i(- + rel)0el(r)C/(r)^o drds 
•) / r e-''(*-^)^^P,(/fi)F(|^^ > M)V2{- - se-{) 

Jt-A Js-B 

0_ei(s)e-^(^-^)^2P,(il2)Fi(- + rel)0ei(r)C/(r)V'o drds 

j'^^T _j_ ^"2, low _j_ j2, high 

Here 5 is a constant independent of Co so that B » A, M » B, and F{\p\ < M),F{\p\ > M) 
denote smooth projections onto frequencies |p| < M, \p\ > M, respectively. The term jf'^ is 
estimated similarly to Jf above. Indeed, using that Vi{- + Tei)0e|(r) = 0ei('7")^i, one has 

II^c'IIl^+L- < \\Jc'^\\L^<f r 1^2(--Sel)0_,l(s)e-(^--)^^Pe(i/2) 

Jt-A Jo 

Vi{- + rel) Qei{'r)U{T)i;o ^ drds 

i-t rs-B 

drds 



< f f 0-ei(s)e-'(^--)^^P,(i/2)0el(r)Fit/(r)V^o 

Jt-A Jo 



< 



< 



t rs-B 1 

||yiC/(T)V'o||LinL2 drds 



t-AJo {s-tYI^ 

\\U{t)iI)o\\l2+l°- drds 



t rs-B 



t-AJo (s-r)"/2' 



^ ^ / * dT ds ,, , ,, /■* dT „ , „ ^ ^ 

^ Co / -. v^/o / ,„/o \m\\L^nL2 + / / 1 rTJ^imhdTds 



t-A JB 



< CoAP-("/2-i)(t)^"/2||V,o||iinL2 +^i3(t)-t||V^ollLinL2. 
Since B >> A and n/2 — 1 > we can conclude that 
(3.22) Pc'IIl^+l- < 10-2Co(t)-"/2||V,o||iini2, 

provided Co is large. 
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3.4 Low velocity estimates 

The idea behind the estimate for Jc ' is that the norm of the operator 

Xi{t, ■)e-<'-'^''^Pc{Hi)F{\p\ < M)V2{- - sei) 
for aU s in the interval [t — A, t] is small. This can be explained as follows: 

The support of xi(t, x) with respect to x belongs to the ball i?2(5t(0)- On the other hand, since V2 
is compactly supported V2(- — se{) is different from on the set Bji{se{). Here R is the size of the 
support of V2, R « to < t. Since s e [t — Ai,t], we have that 

U BrC Br+a{ tei) 

s&[t-Ai,t] 

and thus the support of V2{- — sel) is approximately (1 — 5)t units away from the support of 
Xiit, •)• The operator e~^^^~^^^^ Pc{Hi)F(\p\ < M) can "propagate" the set K into the set L only 

if {t — s)M > dist(i^, L) according to the semiclassical picture. Setting K = sppt^V2(- — sel)^ 
and L = sppt^xi(^) ) we see that this would require that 

{t - s)M > (1 - 5)t. 

However, \t — s\< A and M is fixed so that ^, M << to < 
We make this argument rigorous in the following lemma. 

Lemma 3.3. Let A, M he large positive constants and A, M « t. Then 

sup Wxiit, •)e-^(*-^)^iPe(i^i)i^(b1 < M)V2{- - sel)|U2_^. < 

\t-s\<A Ot 

where 5 is the constant from the definition of xi , X2 • 

Proof. The proof is a commutator argument. Let xi = xit-i^)- Firstly, we claim that 

(3.23) \\[xi,Pcm]h2^L^<e-'^. 

Clearly, 

[xuPciH,]] = [xi,I-Pb{Hi)] = -ixuPbiH,)]. 
Recall that ui, . . . , u-^ are the exponentially decaying eigenvalues of Hi. Therefore, 

ni 

[Xi,-P6(^^'i)]/ = ^ixiMf^n) -Ui{fXi,Ui)) 

1=1 

m 

= ('^"^ + Xi)ui{f, Ui) - Ui{f, (xi - l)ui)) . 

i=l 

In of the support of xi — 1 we have that 

ll(l-XiKI|2<e-^* 
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and thus 

||[xi,n(i^i)]/||L^<e-"^*||/||z.2, 
as desired. Secondly, we claim that 

AM 



(3.24) l|[xi,e-^(*-^)^^]i^(|p1 < M)U.^L^ ^ 



< 



6t 

We write 
and 

Jo 

= ij^ 'e^^^K-VxiV - ^Axi)e-'-^^ dr. 

Observe now that 

|Vxi(t,x)| < ^, |Axi(t,x)| < 

Therefore, 

||[xi,e-(*-)^^]F(|p1<M)||^.^^2 

< It - s| f llVxi Ve-^^'^Fm < M)||i2_^2 + W^xie-^'^'Fm < M)!!^^^^^ 



Since the potential Vi is bounded it is standard that 

(3.25) supllVe— ^VIIl^ < W^fh^ + II/IIl- 

T 

Indeed, 

supllAe-^^^VIlL^ < sup||e-^^^ii/i/||^2 +sup||Fie-^^^VIlL2 < IN^fh^ + II/IIl^, 

r r r 

and 1)3. 25p follows by interpolation with L^. Therefore, 

||Ve-^^^iF(|pl <M)||^2^i2 <M. 
Combining terms we obtain the bound 

llfo..-"-"'-lF(|,1<M)k.^.<^ + ^<2^ 
since M « t, which is (|l-{.24)) . Finally, we invoke one more standard fact, namely 

(3.26) I|[xi,i^(|p1 < M)]h2^L. < M-i llVxilloo < 
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To see this, write F{\p\ < M)f = ['n{£,/M)f{S^)]'^ with some smooth bump function ry. Hence the 
kernel K of [xi,F{\p\ < M)] is 

K{x,y) = M''r]{M{x - y)){xi{x) - Xiiv)), 

and H3.26|) follows from Schur's test. One concludes from estimates 1)3. 23() . 1)3. 24() . 1)3. 26(1 that 

Xi{t, ■)e-'^'-'^'''Pc{Hi)F{\p\ < M)V2{- - sel) 



< 



AM 



It remains to observe that Xi(*; ■)^2(" — sel) = since the supports of Xi(^) ") ^-iid V2(- — sel) are 
disjoint. □ 



We now estimate the term 

rt 



J, 



2, low 



[ r e-'^'-'^'''Pc{Hi)F{\p\ < M)V2{- - sel) g-ei (s)e-^(^-^)^^Pe(i^2) 

Jt~A J s-B 

Vi{- + Te{) Qe-i{T)U{T)4^o drds 
by means of Lemmas 13.31 and 13.21 Indeed, one has 



I t2, low II ^11 t2, low 1 1 



A AT f * f * 

/ 0-el(s)e-^(^--)^^Pc(^2)0el(T)ViC/(T)V'C 



L2 



drds 



^ AM 



\\VlU{T)^l;o\\L2dTds<^\\Vl\ 



A^BM 
^ Vi 



t-A J s-B 

L2nL°= C'oW^^ ||'i/;o||iinL2- 



t-A J S-B 



||[/(t)?/^o||L2+2.°° drds 



Since t >> A,M, 11/^2^2,°° we obtain that 



(3.27) 



2, low 



|l2+l- < W-'Co{t)-''/^ML^nL2- 



3.5 Large velocity estimates via Kato smoothing 

The estimate for Jc ' ^^^^ will require the use of the Kato ^-smoothing estimate. We present a 
variant of it in the following lemma. It differs from the original version (see |CSlj , |KYj , [Sj] , |Vej ) 
in so far as we consider — + V instead of — ^A. Moreover, we employ commutator methods. 
Both the statement and the proof are known, see e.g. |CS2j . |BKj . |Soj . |Doj . so we will only sketch 
the argument. 

Lemma 3.4. Let H = -i A + V , \\V\\oo < oo. Then for all T, R, M > 1, 

l-T 

(3.28) 



sup 

Br 



\\F{\p\ > M)e-^'^fU^BR) dt < C{n, V) ^ 

M2 



L2- 



Here the supremum ranges over all halls Bn of radius R> 1 and C{n, V) is a constant that depends 
only on \\V\\qo and the dimension n. 
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Proof. We first prove the following estimate, which will then imply the lemma. Let tl^{t) = e '^^^ipQ 
with H = -\l\ + V, ||y||oo < oo. Then for all T > and < a, 



' II ^ lloo 
rT 



(3.29) 



sup 



xqgR'^Jo jm" (1 + |x - xoI") 



The multiplier V{V) 2 corresponds to the symbol 2 = ^(1 + 4. It suffices to prove this 
with xq = fixed. The proof is based on taking the commutator of H with m := w{x)x ■ where 



wix) = + 
One has, with ip = 1^(1) for simplicity, 
d 



a> 0. 



dt 



i(mV'(0), V'(O)) - i(m0(r), V(T)) 



([m,y]V',V) + {diiw{x) Xj)-^'ijj,dii;'^ + -{l\{w{x)xj)-^il:,i; 



{-[m,H]m,m)dt 
(3.30) (-[m,//]V,V> 
(3.31) 
(3.32) 
(3.33) 

One now checks easily that the two terms in 1)3. 31() satisfy 

w{x)\{V)~^V'ip{t,x)\^ dx + / {dew){x)xj{\7)'^djij{t,x) {\7y^di'tlj{t,x)dx 

Kv)-lvV'(t,x)|2 



w\V{V)'2^|J\''dx+ / {dew){x)xjdj{V)~2^de{V)~'^ijdx 
- ([(V)^H^V',V(V)-^V.) -([(V)^(5,ii;)(x)x,]^V',5,(V)-^) 
+ -<|A(u;(a;)xj)-^0,v) - (Ky]V',V> 



(3.34) 



> 



w(x) 
1 + 1x1 



|(V)"2Vi/'(t,x)|^dx 



dx. 



Notice that ()3.34|) is precisely the space integral in the desired lower bound from 1)3.29(1 . 

There are several ways to bound the commutators ^1 := [(V)^ ,w] and A2 := [{'V)2,(dew){x)xj]. 
For example, one can use the Kato square root formula as in |JSSj . Alternatively, one can invoke the 
standard composition formula from ^'DO calculus, see |Tay| . This gives Ai = T.^^^ 1 where 

{{\/)^,w} is the Poisson bracket of the symbols (V)2 and w. Moreover, T^^^^i ^ is an associated 

3 

^'DO operator and Ra^ another ^DO operator with symbol ai £ S~2. A similar expression holds 
for A2. One checks that |{(V) 2 , tt;}(^, x)| < |Vi(;(x)| for all x,^. Therefore, the two terms in (|3.32p 
both satisfy 



'A 



V 

(V) 



0(t),v(v)-2V'(t) 



< 



[I + |x|°)-+^ 



+ 



\Ra,V{V)-2mh 



(3.35) 



< c\\mr2 + ^ 



(v)-2vv>(t) 2 

(l + |x|")^+l 2' 
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Finally, the two terms in ()3.33|) are bounded by 

(3.36) (1 + ||F||oo)||^(t)||^ < (1 + ||l^||oo)||V'o||i. 

In the above estimate we have used the boundedness of the multipliers m and A{w{x) xj)-^ on 
L^. Integrating (|3.34() . (|3.35|) . and H3.36() in time, inserting the resulting bounds into (|3.3U|) . and 
finally using that 

\{mi;{0),m)\ + KmV'm, V(T))| < 2||Vo|li, 

one obtains (|3.29|) . To pass from (|3.29jl to (|3.28jl . let XR be a smooth cut-off to the ball Br, so 
that xr tt^-s compact support in a ball of size ~ R^^. Then, by p.29|l with q = 1 

£ \\F{\p\ > M)e-^*^Vo|li2(5^) dt < £ \\xrF{\p\ > M)e-^*^^o|li2 dt 



< 



''\\Fi\p\>M)xRe-''''Mhdt+ r \\[xR,Fi\p\>M)]\\l^,\\e-^'^i^o\\hdt 
Jo 



itHj. ||2 
L2 



< / ||V(V)-5F(|p-1 >M)xiie-*^V'o|li2dt + r(Mi?)-2||V.o||i2 

JO 

< / ||F(|p-1>M)xi?V(V)-5e-*^^o|li2dt + M-i / || [V(V)-i , xr] |lL2l|e-**''V'o 

JQ Jo 

< C{n,V)TM-^R^U4l,. 

The lemma follows. □ 
We then have for Jc ' ^^^"^ 

(3.37) ||j2'i^ist^||i.+^^<||j2>'gh||^^ 

< I r >M)V2(--sel)0_ei(s)e-^(^-^)^^Pe(^2)T4(- + Tel)0el(T)[/(T)V'o||^ 

< r r >M)0_ei(5)V^2e-*(^-^)^^Pc(^2)0ei(r)ViC/(r)^o||L2drds. 

Jt-A Js-B 

We need to commute the cutoff F{\p\ > M) past the first two terms following it. As for the Galilei 
transform, one has 

\\F{\p\ > M)Q,,-^{s)fh = \\F{\p\ > M)e-ie-^fi^/||2 = \\F{\p\ > M)e"V||2 = ||i^(|p-el| > M)fh. 
Therefore, applying this to the final line of (|3.37j) . one obtains 

< f r \\V2F{\p-e^\>M)e''^^-^^'''P,{H2)Qe-Ar)ViUiT)ML^dTds 

Jt-A Js-B 

(3.38) +r r ||[y2,i^(|p-el| >M)] e-^(^-^)^2Pe(^2)0ei(r)Fif/(T)V'o|L2C^Tds. 

Jt-A Js-B 
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The final two terms in (|3.38|) are now easily estimated by means of Lemma 13.41 (integrating in the 
variable s — r) and (|3.26|) . respectively. To apply Lemma 13.41 choose a ball B^i around the origin 
that contains sppt(l/2). The conclusion is that 

Pc' IIl2+L- < \\V2\\l^^ f \\P,{H2)Qei{r)V,U{T)i,^\\L2dT 

Ma Jt-B-A 

r r ||e-^(^-^)^2p^(^2)0ei(r)FiC/(T)V'o|L2dTds 

Jt-A Js-B 

< (miL^^B + M-'B^)\\Vi\\LinL^Co{ty''^^ML^nL^- 
Thus for a sufficiently large constant M one obtains 
(3.39) II J2' i^'gh 11^2+ioo < 10-2Co(t)-"/2||V.o|UinL2. 

Combining (jHUSI), (Tmi) . (1X^1) . (IT^ . (1X771) . (ICTl) yields 



(3.40) 



|Xl(t,-)^'cf/(t)^0||L2+L- < ^W-^/'llV-ollLinL^. 



This is precisely the desired bound on the first channel. Since the second channel can be treated 
by the very same method, see the discussion following 1)3. 10() . only the third channel remains. 



3.6 The third channel 

By (unj, 



+ ^2(- - sel))[/(s)V'o ds 



X3 U{t)tpQ = X3 ^4>Q- iX3 I e 

'0 

where xs = X3(^;^)- By the standard decay estimates for e**"2", 

lle'^^VollL^+Loo < (t)-"/2||Vo||LinL- 



Thus 

IIX3f/(i)^o||L2+L- < (i)""^^ll^/'o||LinL2 + / xse 



L'^+L° 



By Lemma 13.21 and the bootstrap assumption (|2.13|) . 



11(^1 + V2{- - sei))C/(s)^ollLinL2 < (ll^illiinL- + \\V2\\L^nL^)\\U{s)ML-^+L- 

< Co (ll^lllLinLoo + ll^2||LinL-)W""/'llV'0||LinL2- 

Therefore, by the same calculation that lead to 1)3. 14() . 

\\X3Uit)ML^+L^ < ^-"/'lIV'ollLinL^ + Co^-"/'+^(t)-"/'||Vo||LinL2 



+ 



X3e"^'-'^^(,Vi + V2{--sel))U{s)i^ods 



t-A 
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Since A is sufficiently large we only need to estimate the last term, which we can split as follows: 

t-A Jt-A 

It suffices to consider the first term, which is split further by means of the projections Pi)[Hi) and 
Pc{Hi)=Id-Pi,{H{): 



t-A 



t-A 



(3.41) 



+ r X3 e'^'-'^^VMHi)U{s)^Q ds =: £i + £2- 
Jt~A 



Since we are assuming asymptotic orthogonality of U{t)iJjQ to the bound states, see Definition 12.41 
Proposition 13.11 implies that 

ft 



lAI 



L2+LO 



< IIAI 



L2 



< 



t-A 



l^i|lL2nL°°ll-Pb(^^i)t/(s)^o||L2 ds 



(3.42) < ^||yi||iinz.-e-"*/2||V^o||L2 < {ty^/^ML^nL^- 

It is now natural to expand £2 further, cf. (|3.19|) : 



(3.43) C2 = X3 [ e'^'-'^^ViP,{Hi)e-''''^il,ods 

Jt-A 

(3.44) -ix3 f r e'^^-'^^ViPc{Hi)e-'^'-^^"W2{- -Tei)U{T)^odTds. 

Jt-A Jo 

The first integral (|3.43() is controlled by means of the decay estimate for the evolution e~'^^^^ Pc{Hi): 



X3 



ft ^ 

/ e<^-'^^ViPc{Hi)e-''^^%bods 
Jt-A 



, < ri||L2nLoo / ||e--^^Pe(^l)V'0||L2+iood5 



X3 



t rs-B 



(3.45g m\L2nLooA{t)-^^/^ML^nL^ < W-"/'llV'o||LinL2. 

For the second integral (|3.44)) one also invokes the bootstrap assumption (|2.13j) , as well as LemmaES 

e'^^-'^^Vie-'^'-^^"'PciHi)V2{- - Tel)?7(r)V'o drds 
i-t rs-B 

e-ii^-^)iiipjH^)Vo(- - Te;)U(TU<^ drds 

1 



t-A Jo 



<ll^i|lL2nL- r r e-'^'--^''^P,{H^)V2{--Te-[)U{T)i;odTds 
Jt-A Jo 



< 



ll^l||L2nL- 11^211 



'B 



< 



^\\U{t)'iPo\\lHl°° drds 
dr ds 



t-A JO (s-r)"/- 
L2nL- 111^2 IIlIrl- |co J^^^j:^\\ML^nL^ drds 

+ / / 7 — —^ll'^oh^ drds 

Jt-A Jo [S - t)"/^ 
< ||l/l||L2nLoo||y2||LinL-(Ai?-"/2+lCo||Vo||LinL2+AB(t)-t||Vo||L2) 
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since B >> A, and provided Cq is large. As in the proof of the decay estimates for the first channel, 
we split the remaining expression, which is the most difficult part, into low and high momenta: 



[ r e-''^'-'^^Vie-'^'-^'>^'PciHi)V2{- - Tel)[/(T)Vo drds = + ^i^^^^h 

Jt-A Js-B 



X3 

where we define 

= X3 r r e-*(*-^)TF(|^^ <M)yie-^(^-^)^iPe(^i)^2(--re-l)C/(r)V^odTds 

Jt-A Js-B 

C^'^"" = X3 f r e-^(*-^)^F(|^^ >M)yie-^(^-^)^iPe(^i)^2(--re-l)C/(r)V^odTds. 

Jt-A Js-B 

As in Lemma 13.31 one has 

(3.46) sup \\x3e-"^'-'^^Fi\p\ < M)Vih2^L2 < 

\t-s\<A 'J'^ 

This is just a consequence of the following: 

1. sppt(x3(t, •)) n sppt(yi) =0 for t > to = (f )'^" 

2. |Vx3(t,-)|< Ji- 

The first property holds since 

X3{t, x) = 1 - xi{t, - X2{t, x) =0 

on the set 

{x : |x| < 5t} U {x : \x — re{\ < 6t}. 
The bound (|3.46|) leads to the following estimate on C^""" 
II r low 1 1 ^ II /»low|i 

AM /"* 
^ ^^ll^i||L°°ll^2||L2nL°° / / \\U{T)i;o\\L2+L^ drds 

0^ Jt-AJs-B 

a'^mb 

^ — ll^l||L°°ll^2||L2nL°°C'o(i)^"^^||'0o||LinL2- 

Since A, M < B are large but universal constants and t > to = (^)^''", this implies that 

(3.47) ||£ii°-|U2+i.^ < 10-2Co(t)-"/2||Vo||LinL2r(T)V'o||L2+Loo. 

Finally, the high energy part is treated by means of Lemma 13.41 as before. More precisely, one has 

Z:?'^^ = X3 r r e-*(*-^)T[F(H >M),yi]e-*(^-^)^^Pe(^i)y2(--Tel)?7(r)^odTd5 

Jt~A J s~B 

+X3 r r e-'^'-'^^ViF{\p\>M)e-''^'-^^"^P^{Hi)V2{--Te-{)U{T)i^QdTds 

Jt-A Js-B 
_ p high , p high 
— ■'-2,1 "'"''-'2,2 ■ 
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The commutator estimate, see (|c{.26() . 



\\[F{\p\>M),Vi]h2^L^<^\\VVi\\Lo 



implies that 



A R 

(3.48) IIA'f'llL^+Loo < ^Co(t)-"/2||Vo||LinL^ <10-'C^oi-"/'||^o||LinL- 

By the smoothing estimate of Lemma 13.41 appHed to the integration variable u = s — t and with R 
sufficiently large depending on the support of Vi, 



C^'2^''\\l^+loo < I I \\ViF{\p\ > M)e-™^iPe(^i)^2(- - rel)C/(T)Vo||L2 dudr 

Jt-A~B Jo 

^^BR\\V2\\LinLoo f r(r)^o||L2+ioodT< ||yi||L^||y2||LinL-^C'o(t)-"/2||^Q| 

2 Jt-A-B M 2 



M2 

Therefore, with M large, 

(3.49) \\C2'^2^'"\\l^+l^ < lO-^Coll^ollLinL^. 
Combining - we obtain 

which concludes the proof of ()2.14|) , and thus of Theorem 12.71 
3.7 Perturbations of charge transfer models 

It is easy to see that the argument proving Theorem 12 . 71 also yields the following result. 
Theorem 3.5. Consider the equation 

(3.50) -dti; - ^ A^A + V^{x - v^t)^; + Vo{t, x)V' = 

* K=l 

V'|i=o = V'o,2; gM", 

where the charge transfer part is an in Definition \2.1\ and the perturbation satisfies 

sup \\Vo(t, Ollinoo < e. 
t 

Let U{t) denote the propagator of the equation 1)3. 5U() . Then for any initial data -00 £ H L'^ , 
which is asymptotically orthogonal to the bound states of Hj in the sense of Definition \2.4\ ( with 
U{t) instead of U{t) ), one has the decay estimates 

\\U{t)M\L^+L^<{t)-'''^\\M\L^nL^ 

provided < e < eo is sufficiently small with eq independent of ipo. 



LinL2. 
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The proof of this theorem is basically identical with the proof of Theorem 12.71 For example, 
consider again the case of two potentials m = 2. Then the Duhamel formula relative to the first 
channel, see 1)3. 8|) . becomes 

XiPcmumo = Xie-''"'Pc{HMo - iXi f e-*(*-^)^iPe(^i)V'2(- - seI)C/(s)V'o ds 

Jo 

(3.51) -ixi [ e-'^'-'^'''Pc{Hi)Vo{s,-)U{s)i^ods. 

Jo 

Impose the bootstrap assumption (|2.13j) on the evolution U. Then the new term (|3.51j) satisfies 

/■* 3 3 

lldSISIIllb+oo < C{t-s) 2 sup ||Fo(s,-)l|inooCo(s) 2||V'o||in2(is 

JO 0<s<t 

< CeCo{ty^\Min2. 

Modifying the other two channels (|3.9() and in the same way, one concludes that the only 

change from the estimates in the previous section is the addition of a term C e Co{t)~2 \\'ipQ\\i^2- 
Choosing e small, one again arives at the improved bootstrap assumption (|2.14|) for the evolu- 
tion ij{t). 

4 Asymptotic completeness 

Our next goal is to provide the following version of the asymptotic completeness for the charge 
transfer model (see jYa2j . |Gr j ) . In this section we require more regularity of the potential, see 
Lemma 14.21 below. 

Theorem 4.1. Let, as before, ui,... ,Um andwi, ...,we be the eigenf unctions of Hi = — ■y + Vi(x) 
and H2 = ~Y + ^2(2^), respectively, corresponding to the negative eigenvalues Ai,--- , Am o-nd 
Hi, ■ ■ ■ ,fi£. Then for any initial data Tpo £ the solution U{t)ipQ of the charge transfer problem 
()2.9|) can be written in the form 

m £ 

U{t)^o = ^ Ae-'^-V + 5^5,e-^^^*0_e-,(t)u;, + e-*4<^o +7^(^), 

r=l s=l 

for some choice of the constants Ar, Bg and the function (pQ. The remainder term TZ{t) satisfies the 
estimate 

II'^(^)IIl2 — > 0, as t — > cxo 

Proof. We first assume that tpQ E PlL^. The L^-assumption can be removed by L^-continuity of 
the wave operators that appear below. Decompose 

m ■■= = Pb{Hi)U{t)^Po + Pb{H2,t)U{t)i;o + R{t). 

By construction we clearly have 

(4.1) Ph{H2, t)U{t)i;o + R{t) G Ran(P,(Fi)), 

Pb{Hi)U{t)il)o + R{t) e Ran(Pe(^2,t)) 
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We further write 

m 
r=l 

for some choice of unkown functions ar{t). Using ()4.1() we obtain, similar to 

CLr + i (V2(- — tei)^(t), Ur) =0 for all 1 < r < m. 

The exponential localization of u,. implies that |(V2(- — tei)'ip{t),Ur)\ < e^"*. Therefore, ar{t) has 
a limit ar{t) Ar, as i — > +00 and thus. 



(4.2) 



iXrt„ 



Ur 



r=l 



L2 



0, t ^ +00. 



We next define the functions Vr = lim U{t) e * n,.. The existence of is guaranteed by the 

t— >+oo 

existence of the wave operators (see |Grj and for an independent proof. Lemma |4.2|) 

nl = s- limU{t)-^e-'*^'Pb{Hi) 

In this notation, Vr = ilLur- In particular, 



(4.3) 



iXrt„ 



Ur 



r=l 



r=l 



L2 



0, t ^ +00. 



We then infer from (14.21) that 



(4.4) 



r=l 



L2 



0, t^+OO. 



Similar arguments apply to Pi,{H2,t)U{t)ipQ. More precisely, we write 

U{t)i;o = Pb{H2,t)U{t)i;o + T{t) = Q.s,{t)Pb{H2) Qs,{t)U {t)^^ + V(t). 

Therefore, 

(4.5) {t)U{t)^(, = P,{H2) 0el mit)iJo + 0e-i (t)r(t) 

Recall that the function 

^(t) = Q^,{t)U{t)ll^o 

is a solution of the problem 
(4.6) 

According to (ji^J, 7p{t) = Pi,{H2)^{t) +ri(t), where Ti{t) = 0^^(O)r(t). In particular, 

Ti{t) £ Pc{H2)L^ 



-dtip -^iJ + V2ix)i; + Vi{x + tei)^P = 0, i;\t=o = 0e-i(O)V'o 
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Decompose 



s=l 



for some choice of unknown functions bs{t). After substituting the decomposition in (|4.6I) we obtain 
the equations 

bs{t) +i{Vi{- +tei)ijj,Ws) =0 for all 1 < s < £. 

Using exponential localization of Wg we conclude the existence of the limit bs{t) Bg as t ^ +00. 
Thus \\PbiH2)m - ELiBse-'^'^' WsWl^ ^0, t ^ 00. Equivalently, after applying 0-_^i(t), we 
have 



(4.7) 



s=l 



L2 



0. 



We now invoke the existence of the wave operator (see |Grj and Lemma l4.2() 

= s - lim[/(t)^i 0^el(t)e"'*'''n(^2) 

t— >+oo 

which allows us to define ojs ■= u;^. Moreover, 



(4.8) 



U{t) ( Bsu^s) - Yl Bse-'^^' 0_e-i {t)w. 



s=l 



s=l 



L2 



0, t +00. 



It then follows from (|4.7j) that 



(4.9) \\PbiH2,t)Uit)^lJo-Uit){YBsu;s)\\L^ ^0, t ^ +00. 

s=l 

We now define the function 

(4.10) <^ := ^0 - ^ ArVr - Y Bs^s, 

r=l s=l 

which will lead to the initial data (j)Q for the free channel. We have that 

m I 

Pb{Hi)U{t)(t> = Pb{Hi)U{t)ijo - Pb{Hi)U{t){YA^r) - Pb{Hi)U{t){Y,Bsu;s). 



r=l 



It follows from (g^ and the identity P^{Hi) = Pb{Hi) that 



PbiHi)U{t)^o - Pb{Hi)U{t){Y,ArVr) 



r=l 



m I 



L2 



s=l 



0. as t ^ +00 



(4.11) 

Furthermore, 

(4.12) Pb{Hi)YBse''^'^' Q-eAt>j = ^^S,e-^^=*(0_e-,(t)«;,-,ur>n, ^ 



s=l 



r=l s=l 
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in the sense as i — > +00, due to the exponential localization of the eigenfunctions u^. One 
concludes from and (|I?T^ that ||Pf,(i?i)C/(t)0||i2 ^ 0. Similarly, we can show that 

\\Ph{H2,t)U{t)(j)\\i2 — > 0. Thus, U{t)(f) is asymptotically orthogonal to the bound states of Hi and 
H2 and therefore, according to Theorem 12.71 satisfies the estimate 



(4.13) 



\\u{m\L^+L^<{t)- 



In order to be able to apply estimate 1)4. 13() one needs to verify that G L^nL^. By the assumption 
of the theorem the function ^pQ G H L^. Thus it remains to check this property for the functions 
Vr, r = l,..,m and ujs, s = !,..,£. Since Vr = ^lUr and cOs = Ogtfs the property follows 
immediately. The property, on the other hand, is guaranteed by Lemma 14.21 below. Assuming 
this lemma for the moment, we now consider the expression 



U{t)c^ = ^-i 



- {Vi{x) +V2{x- sei)) U{s)(l) ds. 



The estimate 

r°° . A r°° 

/ ||e-*^- {Vi{x) + V2{x- tei))U{s)ct>h^ < mh^nL^ + mh^nL^) / \\U{s 
Jt Jt 

< (||yi||i2nL- + ||^2||L2nL-) ^ 0> as t ^ +00 

allows us to show the existence of the limit 

(j}o := lim e'^^U{t)(j). 

t—tOO 

It follows that 

(4.14) ||C/(t)(/)-e''4^o||L2 ^0, 

Combining (|0|) . and we infer that 



■iXrt ^ 



Ur 



Y,Bse-'^^'5-g,it)ws-e^ 



4, 



L2 



0, as t ^ +00, 



r=l s=l 

as claimed. 

Lemma 4.2. Assume that the potentials Vi{x),V2{x) satisfy the following conditions 



□ 



(4.15) 



E \d2V.{x)\ < 



0<|7|<n+2 



for some positive constants c{k), k = 1,2. Then the range of the wave operators 0,^ ^ ^-^ contained 
in the space of functions. 

Proof. Without loss of generality we only consider the wave operator 

0^ = s - limU{t)-^e-'^^^PbiHi). 

>+oo 
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Therefore, for an arbitrary function / 

m 

^ — * t — >4-nr) 



t—>+oo 
r=l 



where Pb{Hi)f = YlT=i fr'^r for some constants /,.. It follows from the Duhamel formula for the 
equation (|2.9j) generating the evolution operator U{t) that 

U{t)-^e~'*"^Ur =Ur + i [ U{s)-^V2{- - sei)e-''"^Ur ds 

Jo 

(4.16) =Ur + i [ U{syW2{- - sei)e''^'''urds 

Jo 

with (|4.16l) follows since Ur is an eigenfunction of Hi corresponding to an eigenvalue A,-. The 
function is exponentially localized in together with its n + 2 derivatives ^ 

V / e2"l^'l|a>,(x)|2dx < C 

0<|7|<n+2 

for some positive constant a appearing in 1)3. 4() . It follows from 1)4. 15() that the function 

Gr{s,x) := e~*'^''*V2(x — sei)ur{x) 
has the property that for any A; > and multi-index 7, < I7I < n + 2 

\\{x)>'d2Gr{s,-)\\Ll<c{r,h\,k){sr^r^^-'. 

To prove the desired conclusion it would then suffice to show that there exists a positive constant 
k such that for any function g(x) 

(4.17) IKx)[t]+ic/(t)5||^. <(t)3[t]+3 m'd^.ah'^, vt>o. 

|/3|<n+2 

We note here that U{t)g denotes the solution of the problem 

(4.18) -dti^-^il) + V{t,x)^ = Q, V{t,x) = Vi{x) + V2{x-tei) 

evaluated at time t with initial data function g given at time t = 0. The required estimate should, 
however, involve the expression U{t)~^g which stands for the solution of the above problem at time 
t = with initial data g given at time t. Yet it is not difficult to see that the two problems are almost 
equivalent. Therefore, for simplicity we shall prove estimate ()4.17() . We note that the estimates 
of the type H4.17() for problems with time independent potentials are well-known. They have been 
proved in the paper by Hunziker In the time-dependent case the argument is essentially the 
same. More precisely, define the functions 

^.M7|(*) ^= E E \\{^y'd2'U{t)gU2 

j'=0|7'|=0 



^The localization of higher derivatives of Ur follows from the localization of Ur stated in 13.411 and the equation 
-Ur + Vi{x)ur = XrUr with potential Vi{x) which is bounded together with all its derivatives of order < (n + 2). 



28 



for any index j > and any multi-index 7. Using equation H4.18() we obtain that 

Computing the commutator we obtain the recurrence relation 

{x) 



<'^mio) + {tf E m92'v sup <^,_i,|,h.i(t) < 

/l^ol I \ / t-x 0<T<t 

I7 I<2l7l 



^ fc=0 



where C{V) is a constant depending on 



(4.19) 

l7'l<2(|7l+i-l) 



{t) 



T oc 



In addition, differentiating the equation H4.18() (I7I + j) times with respect to x and using the 
standard estimate, we have 

^o,M+j-W < C{V){1 + |r|H+^-)$o,|7|+i(0) 

Therefore, 

$i,H(t)<c(y)(i + |t|)3^-+H<i>,,i^i+,(o). 

Now setting j = [§] + 1 and I7I = we obtain the desired estimate (|4.17|) with A; = [^] + 1. 
Observe that the assumption (|4.15|) controls the constant C{V) in (|4.19|) for the potential V{t, x) = 
Vi{x) + V2{x-tei). □ 

5 L°° estimates 

In this section we develop a general simple scheme allowing to convert the + estimates 
developed above into the true dispersive L°° estimates. 

Proposition 5.1. Let ip he a solution of the Schrddinger equation 

-dti:-^^ + V{t,x)il: = ^, V'|t=o = V'o 
with a time- dependent potential V{t, x) satisfying the condition 

(5.1) sup \\V{t, OIIlihl^ < 00, sup \\V{t, OIIli < 00. 

t t 

Here V{t, •) denotes the Fourier transform ofV only with respect to the spatial variable x. Assume 
that %[) obeys the estimate 

(5.2) llV'COIlL^+ioo <(t)-t||Vo||LinL2. 
Then ^ also satisfies the L°° estimate 

(5.3) \\m\\L^ <\t\~'^\\ML^nL^- 
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Proof. For simplicity consider the case of dimension n = 3. The n > 3 dimensionial case can be 
treated by considering A: = [^] + 2 terms in the Duhamel expansion and repeatedly exploiting the 
cancellation property below (see |JSSj ^. 
We have 



'o 



V{s,-)ij{s)ds 



(5.4) = e 2 ^0 - W e 
/o 



Jo 



t rs 



2V{s, ■)e"^'-^>-V{T, >(t) dT ds. 



We recall the cancellation property 



(5.5) 



sup y-'^V{t, ■)e-''^f\\L. < \\Vit, II/IIlp, yp G [1, oo] 



which was used in |JSSj . This property can be checked directly for pure exponentials, and then 
one writes the potential as superposition of those. We can then easily estimate the first two terms 
in (|5.4|) by the desired bound |i|~2 H'^oIIli- ^ov the second one, split the integration according to 
0<s<l, l<s<i — 1, and t — 1 < s < t. For the values t > 2 we split the term 

/"* r e*(*-^)Ty(s, •)e*(^-")Ty(r, Xr) dr ds = 
Jo Jo 

l-t~2 r-t-l 1-1-2 fT+l i-t-2 ft ft ft 

(5.6) / /+/ /+/ / +// 

^Jo Jt+1 Jo Jt Jo Jt-i Jt-2JT 



y'^('-')TV{s, ■)e'^'-^^^V{T, •)V'(r) 



ds dr. 



With the help of the assumption (|5.2)) we estimate 

ft-2 rt-l ^ ^ 

e<*-«)-T/(s, •)e<"-^)-1/(r, Xr) ds drU'^ 







r+l 



< /■* f' snp, WVjs, -MIli sup, \\V{t, .)hinL- /^^-f ^„^^|L/. n 
Jo Jt {t — s)2 (s - r)2 

(5.7) <{trh\ML^nL^- 

Furthermore, in view of 1)5. 5() . the decay of the free evolution, and 1)5. 2() . 

ft-2 /-T+l ft-2 /■* X A A 

+ / / ]e'^'-'^^V{s,-)e'^'~^^^V{T,-Mr)dsdT 











t-1 



3 



< / 'sup \\V{s, Oll^i ^^P^ (^)-t dr||V^o||LinL2 

10 s 



{t-r): 



< 



{t)-^Uo\ 



LinL2. 



Finally, 



t rt 



t-2 Jt 



5»(*-^)Ty(s, .)e*("-^)Tl^(r, •)V'(r) ds dr 



<{t) 2 sup||y(s,-)||Li SUp||y(T,-)||LinL2 / / \t-T\ 2 ds dr \\lPo\\ ^i^L^ 
s T J t-2 Jt 



(5.8) 



< 



(t) Hl^o||LinL2- 
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The remaining estimate in the case t < 2 can be carried out precisely in the same fashion as in 
()5.8|) with an integral lower limit t — 2 replaced by 0. □ 



6 Estimates for the inhomogeneous equation 

Inspection of the proof of Theorem 12 . 71 indicates that we can extend the L°° estimates also to the 
case of the inhomogeneous equation 

(6.1) -^t^P-^i; + f2vix-v^t)^l^ = F{t,x), ^^14=0 = V'o- 

This will be particularly useful in the forthcoming application |RSSj to the question of stability 
of multi-soliton states. Again, for simplicity consider the case of two potentials, m = 2. Then we 
have the following 

Theorem 6.1. Let be as in Definition \2.1\ and let tp he a solution of equation (|6.1|) with the 
right-hand side F satisfying the condition 

(6.2) |||F||| := sup {^V(T,-)||Li^^r + (t)t+i||F(t,.) 11^2 } <oo. 

Assume that the following decay of the projections of tp onto the bound states corresponding to the 
Hamiltonians Hi and H2: 

(6.3) \\PbiHMt,-)\\L2 + \\Pb{H2,t)^P{t,■)\\L2 < B{t)-^ , 

cf. Definition \2.4\ holds for some positive constant B. Then we have the + L°° estimate 

(6.4) U{t,-)\\L^+L^ < {t)-'^ (WMl^hl^ + infill 



Proof. The proof requires repeating the arguments from Sections |31 As in Theorem 12.71 we in- 
troduce the three channel decomposition and represent the solution via the respective Duhamel 
formulas. The novel terms generated by the inhomogeneous term F can be written in the form 



hit) = XI 



f\-iit-^)iiip^^Hi)Fis,-)ds 
Jo 

hit)=X2Q-eAt) I e-'^'-''>^^P,{H2)Qs,{s)F{sr)ds 
Jo 

l3it)=X3 [ e-*(*-^)TF(s,-)ds, 



which should be compared with (|3.8() . 1)3. 1U() . and (|3.11() . respectively. Observe that by our assump- 
tion (|6.3|) the projections of ip onto the bound states corresponding to Hi and H2 already satisfy 
the desired decay. The terms Ii, 12,13 are estimated in a similar manner. Consider for instance the 
L'^ + L°° norm of /i(t): 

l|/i(t)llL^+L«= <^'(i)-^l|e-'^*-^)^^^'c(^i)F(.,-)||L2+ioo ds + |je-(*-^)^ip,(i/i)F(.,.)||L2d. 
<(t)-t \\F{t,-\\L^nL^ds + \\F{s,-)\\L2ds<{t)-'i\\\F\\\ 
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as desired. 



□ 



7 Decay estimates: Systems with a single matrix potential 
7.1 Three dimensions 

In the following sections we develop the decay estimates for the charge transfer model with ma- 
trix potentials. These systems arise in the study of stability of multi-soliton states on nonlinear 
Schrodinger equations |E,SSj . 

We first investigate the question of the decay estimates for problems with a single time- 
independent matrix potential. We show that the approach developed by Ranch R for scalar 
equations can be naturally extended to systems. The method relies on meromorphic continuation 
of the resolvent of the Hamiltonian across the spectrum and requires exponential decay of the 
potential at infinity. Then, following the ideas we have developed in the previous sections for the 
scalar charge transfer models, we establish the decay estimates for solutions of the matrix charge 
transfer problem satisfying appropriate "asymptotic orthogonality" conditions. Our approach to 
dispersive estimates for systems is more direct than the one of Cuccagna |Cuj . who extended Ya- 
jima's L^'-boundedness result for the wave operators to the case of systems. On the other hand, 
|Cu| does not require exponential decay of the potentials as we do. 

In this section we consider the case of three dimensions, n = 3. The higher dimensional situation 
is discussed in the next section. 

In this section H = —^A + where > and H is an operator on L^(M^). Set 





one can also write 




Since B* 



B it follows that spec(i?) C One checks that for JRz ^ 




where 1)7. 2() also requires the expression in brackets to be invertible, and with 




Vi|2sign(yi) 
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It follows from H7.1|) that spec(-B) = (— ioo, —ifj] U [ifi,ioo) C iR. We will assume that there exists 
eo > so that for all x G 



(7.3) |Fj(x)| < C7e-^«l^l for j = l,2. 

Then it follows from Weyl's theorem, see Theorem XIII. 14 in |E,S4j . and the representation (|7.2)) for 
the resolvent of B, that speCe^g(^) = speCg^g(i?) = {—ioo, —i^j] U [i^, zoo) C iM. This was observed 
by Grillakis [(irrj . Moreover, (|7.2j) implies via the analytic Fredholm alternative that {A — z)~^ is a 
meromorphic function in C \ (—ioo, — i/u] U [i/i,ioo). We need to make several further assumptions 
on A that are motivated by applications to NLS, but we first switch to a different (but equivalent) 
way of writing systems. Let ip he a solution of some NLS. Studying stability questions for ip, as 
we do in RSS^, leads to a system for the variation dip that can be written either for the vector of 
real and imaginary parts of Sip, or for the vector (|^). For the purposes of obtaining (|7.2j) it was 
convenient for us to follow the former convention, but from now on it will be advantageous to use 
the latter. It is clear that these two different ways of writing the system are equivalent by means 
of a change of coordinates in C^. More precisely, one has the following lemma whose proof is left 
to the reader. 

Lemma 7.1. Let 

Then, with H = —^A + fi, 
(7.4) 

holds if and only if 
(7.5) 

where U = l{Vi + V2),W 
Abusing notation, we v 

with U, W real- valued. We now specify the spectral assumptions we impose on A. 

Definition 7.2. Let A = B + V with V exponentially decaying. We call the operator A on 
n := L'^{R^) X L2(r3) admissible provided 

• spec(^) C M and spec(yl) H (— /x, n) = {uJi | < £ < M}, where ujq = and all tOj are distinct 
eigenvalues. There are no eigenvalues in speCggg{A) . 

• For 1 < i < M, Li := ker{A - uif = ker(A - uji>), and ker(yl) C kei{A^) = ker(^3) =: Lq. 
Moreover, these spaces are finite dimensional. 



i \^2j V W -H-U J\ij2 

nil from now on write 

, U -W \ , „ ( H + U -W 
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• The ranges Ran(A - ivg) for 1 < i < M and Ran{A^) are closed. 

• The spaces are spanned by exponentially decreasing functions in TL (say with hound e^*""!^'! ). 

• All these assumptions hold as well for the adjoint A*. We denote the corresponding (gener- 
alized) eigenspaces by L|. 

• The points zizfi are not resonances of A (we will define a resonance below, see Remark \7.1U<{j . 

Definition 17.21 is motivated by applications to the questions of stability of soliton solutions of 
NLS (see e.g. (RPT] . (ESS])- 

Henceforth it will be automatically assumed that A is admissible. 
Lemma 7.3. There is a direct sum decomposition 

M M 



(7.7) w = E^^ + (E^I 



3=0 j=o 

This means that the individual summands are linearly independent. The decomposition (|7.7|) is 
invariant under A. Let Pc denote the projection onto ^X^jlo-^j) which is induced by the split- 
ting 1)7. 7|) . and set = Id — Pc. Then APc = PcA, and there exist numbers Cij so that 

(7.8) Pb/ = ^0,Q,•(/,V'^) for all f^H 

where 4>j,il^i are exponentially decreasing functions. 
Proof. One has 

M j_ M 

(7.9) (^Yl ^) = Ran(A2) n f] Ran(A - ujj) 

j=0 j=l 

by the assumption that the ranges are closed. We verify first that the summands in ()7.7|1 are 
linearly independent. To this end suppose 

M M 

(7.10) / = ^ /,-, eL„0<j<M, / G ( i*) . 

j=0 j=0 

Then by ((7^ . / = A'^go and f = {A- u)j)gj for some go G Dom.{A^), gj £ Dom(A). Note that for 
j >l, fj ker(^ -LUj) ^ fj e Bom{A^) for all s > 1. Thus (TTTT]!! imphes that 

M M 

A^f = Y,A^fj, and A^ll{A-ujj)f = 0. 

i=i j=i 

Therefore, by our assumption on (generalized) eigenspaces 

{A - ioef lliA - iOj)A^ge = ^ {A - u;i) JJiA - ujj)A^g^ = 0, 
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and similarly, 

M M 

A'' 1[{A - oj,)go = O^A^ll{A- uj,)go = 0. 
Hence / satisfies the equations 

M 

Yl{A-LOj)f = 0, A'^Yl{A-ujj)f = for all 1<£<M. 
3=1 j¥« 

Using that Lj n = for j ^ k, one concludes inductively that nj>jo(^ ~ ^j)f ~ *-* Jo- 
Indeed, setting Aq := and Aj := A — uj, we just showed that Ylj^iAjf = for all < £ < M. 
Thus also Ylj^£ k^jf ~ ^ i ^ k. Continuing inductively one finds that / = 0. By the 

linear independence of the spaces {LjI^Iq this implies that /o = /i = • • • = Im = 0. Thus, 

M M ^ 

j=0 i=0 



(E 



Dually, one finds that 

M , M 



(E 



j=0 j=0 

which means that the right-hand side of ()7.7() is dense in 7i. Since the sum of a closed space 
with a finite dimensional space is again closed, we have proved that (|7.7|) is indeed a direct sum 
decomposition. Hence every f ^ TL can be written as f = g + h where g G Ylj=o^j^ ^^'^ ^ ^ 



EjLoLj) ■ We set PJ = h, P^f = g. Since ALj C L, for < j < M, and A((Ei=o^jj ) ^ 

M T *\ /ID DA 4-„ — 4-^ ,„„„! :„„4-: — — j-„ v~^M 



Sj=o -^i j ' '^^^ ^-^c = ^c^- Denote the orthogonal projection onto Ylj=o Q- Then 

there is an isomorphism 

Qf ^ Pbf 



S : 



i=o ~^ Z^j=o 



and Pfy = SQ by construction. As a map between finite dimensional spaces 5* is given by a 
matrix {cij}. By assumption, these finite dimensional subspaces of W are spanned by exponentially 
decaying functions, and H7.8() follows. □ 

Observe that H7.2|) implies that there exists A > large so that 

\\{A- {iX + z))-'^\\ <\Qz\-'^ for all Qz>0. 

Unless indicated otherwise, || • || refers to the operator norm on 7i. This allows one to define the 
(quasi-bounded) semigroup e**"^ by means of the Hille-Yoshida theorem and He**"^!! < e*^. We now 
make one more requirement for A to be admissible, namely the stability assumption 

(7.11) sup ||e**^ Pell < oo. 
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Under this assumption the semigroup is bounded on all functions that have no component in the 
space Lq under the splitting (|7.7() . Otherwise, it can grow at most like t. 

Assumption ()7.11|) is related to the notion of linear stability in the context of stability of soliton 

solutions of NLS (see e.g. |Weinj . ^^). 

Our goal in this section is to prove the following theorem. 



Theorem 7.4. 

• Let A be admissible as in Definition \7.^ together with the stability assumption 1)7. llf) . Then 

(7.12) \\e''^Pcf\\L^+Lo^<{t)-h\f\\L^nL^ 
for all t. 

• If in addition the matrix potential satisfies ||V^||li < c>o, then one has the bound 

(7.13) l|e**^Pc/||L-<ir^||/||LinL^ 
for all t. 

The norms here are given by Definition 12.31 

:= inf/=h+s(||/i||L2 + ||flf||L°°) 

with the dual space (L^ + L°°)* = n L^. 

The proof of Theorem 17.41 will follow an old strategy of Ranch ^ , see also jDMcLT] and |Vaij , 
that is based on representing e**^ as the inverse Laplace transform of the resolvent, and then 
to move the contour across the spectrum. This only gives local decay, and we then use an 
observation of Ginibre [(irinj to pass from this to the bound (|7.12|) . In the following we prepare the 
proof of Theorem 17.41 by means of several technical lemmas. We will assume from now on that A 
is as in Theorem 17.41 

The next lemma is a version of von Neumann's ergodic theorem. 



Lemma 7.5. For any u) G 



1 



(7.14) ^ I e-'^^h'^^Pcdt ^0 asT^oo 

T Jo 

in the strong sense. 

Proof. Let {ujj}j'LQ be as in Definition 17.21 We first prove 1)7. 14() for lo equal to one of the u>j, j ^ 0. 
Let / G RanPc- Then for every 1 < j < M there is gj £ Dom(^) so that f = (A — ujj)gj, see (|7.9|) . 
Note that gj is unique up to addition of an element of Lj . Thus 

(7.15) ^ £ e— .*e^*^/ dt = ^ (e^^^^-^)^^ - 9j) ■ 
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It remains to show that we can choose gj so that H7.11() apphes. To do so, use H7.7() to write 
gj = Ylk=o + h where G and h G (X^fclo -^fc)"*" ~- ^- assume that hj = 0. Then 

Agj = Ylk=o^^k + Ah and by A-invariance of (|7.7|) . Ahk G Lk, Ah £ L. On the other hand, 



M 



Agj = f + ujjgj = ^ ujj hk + ujjh + f. 



k=0 

The last two terms lie in Ran(Pc) = L, whereas clearly ujj £ for all < A; < M. By uniqueness 
of such a representation, Ah^ = ujjhk for all < /c < M. But then 

hk G ker(^ - Wj) n ker(A - LJfc) = {0} for all k ^ j, 
ho G ker(A - ujj) D ker(^^) = {0}. 

Since hj = by assumption, it follows that gj = h G Ran(Pc)) as desired. Hence the right-hand 
side of ()7.15|) is 0{T~^) and we are done with this case. Now consider the case to = wq = 0. We 
write / = A^g for some g G Dom(74^). By (|7.7|) one has as before Ag = Xlfclo ^k + h where hk G Lk 
and h £ L. We claim that — + Ag G Ran(Pc) and Aho = 0. If so, then one has / = A{—ho + Ag) 
which implies that 

1 r'^ 1 



T 



^ e^'^f dt = ^ [e'^^{-ho + Ag) - {-ho + Ag)) , 



cf. (|7.15j) . and the right-hand side is again 0{T ^) by (|7.11|) . To prove the claim, simply observe 
that 

M M 

f = A'g = Y,Ahk+Ah^Yl G L. 

k=0 k=0 

Since Ahk G Lk one concludes from Lemma 17.31 that Ahk = for all < A; < M. If 1 < A; < M, 
then hk G Lq D Lk = {0} so that hk = 0. But this is precisely the claim, and we are done with the 
case w = as well. 

It remains to consider u; G M \ {wjljlg- Firstly, if a; G (— /i, /x), then {A — u!)~^ is a bounded 
operator on 7i by assumption. Hence, 

i £ e-*-*e^*^Pe/ dt = {A- u;)-'l,(^e'nA-^)pj _ pj 

and the right-hand side is 0{T^^) by (|7.11jl . Finally, assume that |ti;| > /Li so that lo belongs to the 
essential spectrum of A. This case is different from the previous ones because one cannot expect 
to obtain the 0{T^^) bound. Since 



sup 

T 



T 

^-iu^t^itApj 



< oo 



for all / by (|7.11j) it suffices to prove that 



1 

T 



T 
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for all / belonging to a dense subspace of Ran(Pc). In view of the preceding, it therefore suffices 
to show that for all > /u, 



(7.16) 



{A - a;)Ran(Pc) = Ran(Pc 



By the A-invariance of Ran(Pc) it is clear that (A — a;)Ran(i-*c) C Ran(Pc). By Lemma l7.3l we can 
write a direct sum decomposition Ti = L' + L where L' = X^jlo ^"^^ ^ space from above 

with finite codimension. Now ker(A* + ijo) = {0} implies that 



(7.17) 



n = {A- uj)n = {A- uj)L + {A- uj)L'. 



Note that by j4-invariance of the splitting both summands here lie again in the spaces L and L', 
respectively. In particular, the sum on the right-hand side is again direct. Since ker(^ — ijo) = {0}, 
one has dim((A — lo)L') = dim(L') < oo so that {A — lo)L' = L' . By the directness of the sum 
in (|7.17|) therefore 

{A-uj)L = L, 

as desired. 



□ 



The last of the abstract statements is a bound on the resolvents. 



with > one has 



Lemma 7.6. For all z S 

(7.18) \\iiA-zy'Pc\\< 

for some constant C. 
Proof. Since for "Siz > 0, 

poo 

{iA - z)~^P, = - / 
Jo 

()7.18|) is an immediate consequence of ()7.11() with C 



C_ 

Wz' 



suPigM||e**^Pc| 



□ 



Following Ranch j^, we now show how to continue the resolvents {iA — z)^^ = —i{A + iz)~^ 
meromorphically across the spectrum as bounded operators in an exponentially weighted L?' space. 
Set {Ei;f){x) := e~^P^^^ f {x) where p is smooth, p(x) = |x| for large x and e > will be some small 
number (compared to eo in l|7.3j) above). Our goal now is to prove that E^liA — z)~^Ee has a 
meromorphic continuation to the region which lies to the right of the curve shown below. From 
now on, we set /i = 1, so that H = — + 1. Moreover, we denote the curve depicted in the 
following figure by Fg and the open region to the right of it by T^r^ • 




/\ 
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We start with the resolvent of B. By H7.6() one has 
The resolvents 

I p-^2(l±iz)\x~y\ 

(7.20) {H±iz)-\x,y) = {--A + l±iz)-\x,y) 



2 2Tr\x — y\ 

have singularities at the points ±i, respectively. We will resolve these singularities by means of the 
transformations Ui : C, ^ z = i — iC,"^, and U2 '■ C, ^ —i + iC,"^- For the domain of C/i we will choose 
either of the regions 



(7.21) 7^l=7^l(e) := {cGC|3f?C>-^, ^C>-^, -^<ArgC<^, |Cl<r} 

(7.22) 7^2 = 7^2(£) := {c G C | 5RC > J < Arg C < 



Here e > is a small parameter, and r will be chosen to be on the order of 1. These two regions 
are shown in the following figures. Firstly region TZi, with on the left and z = C/i(C) on the right, 



and secondly region with Q on the left and z = Ui{C,) on the right: 




The curved (but not circular) pieces of the boundaries that are shown in the z-plane are parabolic 
arcs which can be of course written down exactly. For the map U2 we choose the domains 1Z\ and 
7^2 which are the reflections of T^i and respectively, across the real axis (i.e., conjugation; note 
that U2{0 = t^i(C))- now defines T^r^ H {z G C 1 5Rz < 1/2} by means of four separate regions, 
namely by ?7i(7^i), C/i(7^2)i t^2('^i)) ^^2('^2) with r sufficiently large, say r = \/5/2 (by taking 
£ small we can still ensure that we move only a very short distance into the negative half plane 
< 0, at least on compact sets). This choice of r is made to ensure that the two disks intersect 
the real axis without coming close to each others centers. 
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Lemma 7.7. Let < e be small. Then the weighted resolvent E^iiB — z) ^E^, originally defined on 
?R.z > 0, admits an analytic continuation to the region with values in the compact operators 

onTL. This continuation, which we denote by F^(z), satisfies the estimate \\F£{z)\\ < Ce(l + |2;|)~2 
for all z £ TZr^^^. Furthermore, Ffr{z) is analytic on the union of the regions z = i — iQ"^ and 
z = -i + iC^ for Id < |. 

Proof. It was shown by Rauch that (recah that H = —^A + 1) 

c^E,{H+ i(i - iOr^E, = E,{-^A + er'Ee 

has an analytic extension to the domain 3^^ > — e/4 with a norm bound of (1 + |C|)~^ (see equa- 
tion (2.12) in [E|). Similarly for ( ^ Es{H - i{-i + iC^))-^Ee = Es{-^A + C^y^Ee. This shows 
that both E^i^H ± iz)~^E^ have analytic continuations to the region 'R-r^^^ satisfying the norm 

bound (1 + as \z\ — > cxd in TZr^^^- The lemma follows. □ 

We can now formulate our meromorphic continuation result for the perturbed resolvent. 

Corollary 7.8. Let < e be small depending on Eq from Definition \ 7. ^ Then for some S > 0, 
the weighted resolvents E^{iA — z)~^PcE^, originally defined on SRz > 0, admit a meromorphic 
continuation to the region TZr^ with values in the compact operators onTi. This continuation, which 
we denote by Ge{z), has only finitely many poles in TZr^ all of which lie in 5Rz < 0. Moreover, it 
has continuous limits on the two horizontal parts ofTs \ {±^} from above and below. Finally, for 
large \z\ there are no poles and one has the estimate 

(7.23) \\Ge{z)\\< 



for all large z € TZ-pg ■ 

Proof. Since Pc = Ld — Ph and one has the explicit representation 1)7. 8() with exponentially decay- 
ing it suffices to prove the same statement for E£{iA — z)~'^Ei, = —iEi,{A + iz)~^E£. By the 
resolvent identity, 

Ee{A + iz)-^Ee = E,{B + iz)-^E,-E,{B + iz)-^E,E-^VE-^EM + iz)-^E, 
Ee{A + iz)-^Es = (l + EsiB + izy^EsE-^v) ^ Ee{B + iz)-^ E^. 



Lemma 17.71 and the analytic Predholm alternative imply that E^^A + iz) ^E^ has a meromorphic 
extension to T^r^/^ which is analytic if \z\ is large. Moreover, 



(7.24) G,{z) = {1 + iF,{z) E-'Vj F,{z) 

for all z £ T^r^^^ at which the inverse is defined. Hence (|7.23|) follows from the previous lemma. 
The poles can only accumulate on the boundary of Tlr^/^- If (5 < |, then the poles of Gg can 
therefore only accumulate on the two horizontal pieces of Ts- We claim that reducing 6 even further, 
one can eliminate any accumulation on these horizontal pieces as well. To see this, "resolve" the 
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singularities at the points ±i by means of the maps z = Ui{C) = i — iC^ and z = U2{C) = —i + iC'^, 
respectively. Clearly, 



(7.25) G,(C/i(C)) = 1 + F.(t/i(C)) E-'V) F,(C/i(C)) 



as long as f/i(C) C TZr^^^- In particular, this holds for all C £ 7^(e/4). In view of the final statement 
of Lemma 17. 7| one can apply the analytic Fredholm alternative to (|7.25|) with |C| < |. Thus 
Gs{Ui{C)) is a meromorphic function for |C| < |. It follows that for (5 > small, G£(f7i(C)) is 
analytic on |("| < (5 with the possible exception of a pole at C = 0. In the z variable this means 
that Gsiz) has continuous limits (from above and below) on the horizontal piece of Ts emanating 
from +i, with a possible pole at z = i. The same statements also hold at —i. □ 

We are now ready to prove the local L^-decay. 

Proposition 7.9. Let A be as in Definition \ 7. ^ together with H7.11|) . Then for e > small one 
has 

(7.26) \\Ej'^P,Ee\\<{t)-l 
for all t. 

Proof. Suppose t > 0. By the inversion formula for the Laplace transform (to be justified below), 

E,{e''^ - I)P,E, = — / e'^iESA - t)-^Pc E, + E,P,E,t-^) dr 

for any a > 0. The subtraction of the identity serves to make the integral absolutely convergent. 
Indeed, {iA — T^{iA — r)^^ = Id implies that 

E,{iA - t)-^PcE, + -EePcEe = -EeiiA - T)-^P,E,EjHAEe. 

T T 

By means of a commutator calculation one checks that E^'^AE^ is a bounded operator W'^''^ x 
W'^''^ 7i. Hence, by H7.23() . the left-hand side decays like |t|~2. We now deform the contour to 
as in Corollary EHl This gives 



Ee{e''^ - I)PcEe = 1^.1 e'^ (Ge{T) + EePcEeT-^) dT 



- J2 ^esi^e'^(G,{T) + E,P,E,T~^y,T = Tj 

Tj poles 

This can be simplified to 

(7.27) E,e''^P,E, = ^ ^ e*" G,(r) dr - Res(e*^G,(r); r = r,) . 

Tj polos 

The integral over converges by (|7.23() . Clearly, 

^ Res(e*^G,(T);r = T,) = ^ p,(^)e*^^ 

Tj poles Tj poles 
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where Pj{t) is a polynomial with coefficients given by compact operators of degree equal to the 
order of the pole minus one. As far as the integral over Ts is concerned, only the horizontal 
"thermometers" can contribute much. "Thermometer" here refers to the fact that on the horizontal 
pieces one first moves to the right on the lower edge, then around a loop encircling iti, and then 
to the left on the upper edge. It is clear that the parabolic arcs give an exponentially small 
contribution. Denote the two thermometers around i and —i by Ti and T2, respectively. To 
determine their contributions, we need to again use the (^-variables. Indeed, setting z = Ui{C) = 
i — iC'^ we showed above that is analytic for z = i — and small ( up to possibly a pole 

at 0. Writing the Laurent expansion Ge(i — i^^) = YIJL-n -^jC"'' one deduces from Lemma 17^ that 
N <2. Also, our assumption on the absence of resonances means that = 0. Therefore, 



(7.28) 



1 
27ri 



I J J- 2m Jj- T -I 27n Jj- 



where ||C(r)|| < |r — i|2 so that the contribution of the second integral is 0(|t| 2). The first is 



equal to the compact operator e 
estimate 



The conclusion is that (|7.27|) now takes the form of the 



(7.29) 



Eee 



itA 



j 



< \t\ 



where we have only retained the purely imaginary poles Tj = e^^^ with distinct rjj. In view of (|7.11j) 
it follows that for large i. 



sup ^Pj(i)e* 



< 1 



which implies that necessarily each pj = const. Finally, the ergodic average lemma. Lemma 17.51 
shows that pj = 0. Hence we are left with 



<|t|-5. 



for large t as desired. 



□ 



Remark 7.10. According to Corollarv 17.81 fand its proof), Gs{i — iC'^) and Ge{—i + iC'^) can at most 
have poles of second order at C = 0. It is well-known that in generic situation there will be no 
singularity at the origin. In fact, a second order pole means that ±1 are eigenvalues, whereas a first 
order pole is what one calls a resonance. For more on this, see for example |R] . 

We now pass from local L^-decay to global decay. 

Proof of Theorem \7.4\ We start with the first statement in the theorem. By Duhamel's formula, 
and with A, B, V as in (|71t|) . 



(7.30) e'^^P, 



ft ft 

, , (,i{t-s)B y^isAp^ ^ ^itBp^ _ ■ / ^i{t-s)B yp^^isA 

lo Jo 



(7.31) Pce''^ = Pce''^ + i f Pj'^'-^^^ Ve'^^ dr = P^e''^ + i f e^^^"")^ PcVe'^^ dr 

Jo Jo 
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Jo 



At-s)B y^i{s-r)A p y^irB ^^^^^ 



Inserting ((73T|) into (f73n|) yields 

ft 

^i(t-s)B yp^^^sB ^ 

Therefore, applying Proposition 17.91 to the middle term in the double integral, yields 

+ / (i^-'S) 2(s) 2 ||^-Pc||L2+L°°^LinL2 IIV'ollLinia 
JO 



+ 



{t-s) 2{s-t) 2{t) 2 drds \\E^^V\\l2^L^nL2 imilL^nL^ 



'0 JQ 



(7.32) 

To pass to 1)7. 32p one uses that 

\\yPc\\L'2+L°°^L^nL2 < ll-Pc||LinL2^LinL2 < OO, \\VE~-^\\L2^Lif^L2 < OO, 

which follow from 1)7. 8p and the assumption that V,(l)j,ipi are exponentially decaying. 

To remove we use the same device from |JSSj that already appeared in Section [3 for the same 

purpose. However, some care is needed for the case of systems since the analogue of 1)5. 5() . viz. 



sup We^^^Ve '"■""llp^p < OO 

l<p<oo 



itB\ 



does not hold now. Indeed, by 1)7. 6() . and with H = — + /i, 

/ i^itHjjp—itH 

(7.33) 



and the off-diagonal terms are not bounded on L*' if p 7^ 2, W ^ 0. As in (|5.4j) . 

JtA2 JtB„ 



tHjj^itH ] 1 



(7.34) 



its T 







Jo Jo Jo 



Jit~s)By^i{s-r)By^trA^^ 



In what follows it suffices to assume U = 0, since the diagonal terms in l|7.33jl have the same 
cancellation as in the scalar case, see 1)5. 5p . Consider one of the off-diagonal entries of the second 
term in 1)7. 34() for times t > ^. Thus, 



e'^'-''^^We-''^^fods 



< 



+ 



+ 



1: 



.1 

< l'\t- 2s\-l lle^^^I^e-^^^lli^illVi'olli ds + 

ft 



\t — s\ 2 \\W\\i \s\ 2 ||'0o||i ds 



(7.35) < (||l^||i + ||VF||i)|t| 



^o||i, 
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where we have employed the cancehation property H5.5|) twice. For smah times < t one obtains 
similarly 



< 


I 


oo 





+ 



< 



\t - 2s|~3 ||e'^^VFe-^"-^||i^i||V'o||i ds + 



|-S|<t2 



\t-s\-l \\W\\i\s\-l\\i^o\\ids 



+ 



(7.36) <m\\i + \\w\\im-'\mi- 



We now deal with the third term in (|7.34|) . and we select one representative off-diagonal term of 
the form 

rt fs 



(7.37) 




JO 



where iPi{t) is the first component of the solution 'iP{t) := e'^'^^ipQ. For t > 4, introducing the new 
variable r' := 2s — r, we rewrite ()7.37() in the form 



/ / e'^'-''^^We'^'-^"^"U^i {2s - t') dr'ds = 

Jo Js 

rt—1 p2s rt p2t—s ft p2s \ 

/ / +/ / +/ / )e*(*-^)^W(^-"')^C/V'i(2s-T')dT'ds 

Jo Js Jt-lJs Jt-lJlt-sJ 



(7.38) 

The estimate for the first term in ()7.38|) is straightforward 











{t-sY 



ir\2 



in2 



(2s-r')-^ll^o||in2drds< 



(7.39) 11^^011102 

The second term in ()7.38() is dealt with in the following manner: 



t f2t-s 



t-1 Js 



e'{t-s)H^^^i{s-r')H^^_^^2s - r') dr'dsWoo < 

ft f2t-s 



t~l Js 



\t - s\ 



< 



(7.40) 



{t)-l 



(t) 2||^l|2||^||2noo 

t 



2||i^ ||2noo(2s - r') ■2\\^po\\ln2 dr'ds 

t r2t-s g 

\t - s\~2\\lPQ\\if^2dT'ds < 



t-1 Js 



1 — » 3 

\t - s|"2 11-0011102 ds < t"2 



yo||in2 



t-1 
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The last term in ()7.38|) requires further sphtting 

ft f2s 



/ / e'^*-'^^We'^'-^">"UM^s - t') dr'ds = 

Jt-l Jlt-S 

Jt-l J2t-s Jt-1 Jt+^ J 

Observe that in the region (s, r') G [t — 1, t] x [2t — s, i + |] we have 2s — t' > On the other 
hand, in the region {s, t') G [t — \,t\ x [t + |, 2s] we have |t — t'| > | — |. Therefore, 



t-l J2t~s 

t rt+ 



I I 2 ||e^(*-«)-^PFe-^(*-^)-^||oo-.oo||e*(*-"')^C/'0i(2s - r') dr'ds||oo dr'ds < 
f \t-T'\-'2{2s-T')-HM\m2dT'ds< 

J t-l J2t~s 



i||t^llin2 



f r"\t-r'\-l\\Min2dT'ds< 

Jt-l J2t-s 

(7.42) , (li-^r^ -|5/2r^)||V^o||in2ds<ri||V^ol|in2 



Similarly, 



e'(*-^)^W(^-^')^C/Vi(2s - r') dr'dsHoo < 

i-l Ji+f 

t r2s g g 

\t - t'\—2{2s - r')-2 ||V^o||in2 dr'ds < 

2 



i||t^llin2 



t-i Jt 



rt r2s 

(7.43) ri / / (2s-r')-^||V^o||in2dT'ds <t-i||V;o||in2 

Jt-i Jt+^ 

For the values of t < 4 we decompose (|7.37|) . with the variable r' = 2s — r, as follows. 

»i r2s 



e 



'''(t-')Hwe'^'-^"^"UM'^s - t') dr'ds 



Js 

.3t 



4 4 
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Note that {2t — s) < 2s provided that s > In particular, this holds for s > ^. We have 

Jo Js 



2s 



\t-s\—^\\W\\2\\Uh{2s-r')—2Uo\\in2dT'ds 



<t--2\\W\\2\\U\\2 r / {2s-T'rlm\in2dT'ds< 

Jo Js 

(7.45) t-IUo\\in2 

Now observe that in the region (s,r') G [f ,t] X [s,2t - s] the value (2s - r') > |. Therefore, for 
the second term in (|7.44j) we obtain that 

ft r2t-s 

(■t f2t-s 

\t - 51-2 ||VF||2||[/||2noo(2s - r')-2 ||V^o||in2 dT'ds 

ft f2t-s 

3 ,, -> 





< 




{t)-2\\W\\2\\U\\2 / / \t-s\-2m\in2dT'ds< 



f 3t 
4 

3 / * , , 1 ,, r ,, , ^ 3 ,, 



(7.46) / |t-s|-2||V'o||in2(is <t-2 



yo||in2 



Similarly to H7.41() we split the last term in H7.44() 

ft r2s 

gi(t-s)Hp^gi(s-r')/^f;^^(2s - r') dr'ds 




2t-s 

(pl pt-\-— pt c2s \ 
/ / '+/ / )e*(*-")-^W(^-^')-^C/V'i(2s-r') dr'ds 

Observe that in the region (s, r') G t] x [2f — s, t + |] we have 2s — r' > |. On the other hand, 
in the region (s, r') G [^, t] x [t + |, 2s] we have |i - r'| > ^. Therefore, 

' ei(*-«)^^VFe^(^-"')^?7V^i (2s -r') dr'ds I loo = 

J2t-s 

"'^^ ||e^(*-^)^Tye-*(*-^)^||oo^oo||e*(*-"'^''f/V'i(2s-r')|U dr'ds < 

J2t-s 





r 3t 



i||f^llin2 




|t - r'|-2 (2s - r')-2 llvJollins dr'ds < 



2t- 



3 /■*+-2, „ 3 



(t)-2 / / |t-r'|-2||V;o||in2 dr'ds < 



(t)"^ fl^-^r^ - |s/2rM ll^o||in2 ds <t-t||Vi'o||in2 
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Similarly, 



Il^lli||t^llin2 / / \t-T'\-l{2s-T')-IUo\\in2dT'ds< 



3 r^' „ 3 „ ^ 3 , 



at 



t'-2 / / {2s-T')—2\\^4ir^2dT'ds<t~-2\\^o\\in2 



The conclusion from the preceding is that 

(7.48) ||e**^PeV^o||oo = \\Pce''^Moc < l^l" ^ ||V^o||in2. 

as desired. □ 
Remark 7.11. Dualizing (|7.48j) yields 

||e**^>,>o||2+oo<ir^llVi'o||i. 



By construction, Ran(P^*) = (kerPc)^ = (ELi^j) > ker(P*) = (RanPc)^ = ELi ^^ and 



{P*) = P*- It follows from these properties that P* is precisely the projection onto (^Ej=i 
corresponding to the decomposition (|7.7|) with the roles of Lj and L*^ interchanged. Therefore, one 
has the corresponding estimate for e**"^ as well, namely 

||e**^Pe^0||2+oo < \t\-'^\m\l- 

Remark 7.12. We have made no use of the wave operators Vl := limt^oo e^^'^^e**^ in this section. 
It is easy to see from the dispersive bound for the free evolution e**^ and Cook's method that this 
limit exists. There is also a corresponding asymptotic completeness statement, namely that is an 
isomorphism from 7i onto RanPc. A first step in this direction is to check the existence of the limit 
limt^oo e~**^e**'^Pc. This is nontrivial, and can be obtained by Cook's method using Theorem l7.4l 
Then there is an additional issue of showing that Ran(r2) = Ran(Pc). For this and much more on 
the wave operators see |Cuj . Let us just mention that Lemmas 17.51 and 17.61 are both immediate 
consequences of this asymptotic completeness. Indeed, Lemma 17.61 follows from 

(A - z)-^p, = n{B - zy\p,ny^p„ 

whereas Lemma 17.51 is an immediate consequence of the fact that 

^ = £ e-'^^'e'^iPcny^Pcdt 

via the usual ergodic theorem for the unitary group e*^. But since the proof of asymptotic com- 
pleteness requires Theorem 17.41 we needed to give a direct proof of these lemmas. 
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7.2 The case of higher dimensions 

The proof that was just presented for n = 3 also works in odd dimensions n > 3. We do not address 
the situation in even dimensions here since it is not clear at the moment whether these methods 
can be extended to that case. 

There are three places where the dimension becomes relevant in the previous proof. The first 
instance is Lemma 17.71 which establishes the analytic continuation of the weighted free resolvent 
E£{iB — z)~^E£ from JRz > to the region TZre , together with the bound (1 + |2;|)~2. Since 

+ ^^-[ E,aA-l + izr'E, 



this reduces to the same statements for the scalar operators (— + 1 + iz)~^ = (— + C^)^^ 
and (^A — 1 + iz)~^ = (— |A — C^)~^, where we have used the maps z = f7i(C) = i — iC^, 
and z = U2{C) = —i + iC'^, respectively. As far as the region is concerned, observe that TZr^ C 

C/i(7^i U7^2) nC/2('^i U7^2)- Rauch proves the desired bound on the analytic continuation of the free 
resolvent by means of the sharp Huyghens principle in dimension n = 3. The very same argument 
also applies to any odd dimension n > 3. Indeed, one writes^ 



(7.49) 



;>oo 

/ e""^* Ee cos{t^/^)E, dt = CE,{-A + C^)"^^£, > 0. 

Jo 



If n is odd, then the sharp Huyghens principle implies that cos{ty/—A)Ei;\\2-^2 ^ e^^*/^ for 
large t, so that 1)7. 49(1 allows one to analytically continue the right-hand side to 5RC > — § together 
with the desired bound. See ||R| for details. 

The second place where the dimension becomes relevant is in the proof of local decay, see the term 
involving C{t) in (|7.28j) . In general (odd) dimensions, there is still the Laurent expansion 

oo 

G,{z)= ^ B,{z-i)i, 

and similarly about the point z = —i. By our assumption on absence of eigenvalue and resonance 
at i one has B-2 = B-i = 0. In order to obtain the desired decay estimate on the evolution e**"^ 
one needs to show that the coefficients Bj of odd powers j vanish for 1 < j < n — 2. This follows 
from the corresponding property of the "free" function F^{z). Indeed, by the previous paragraph, 

°o 2 

(7.50) Fe{z) = Y,Cj{z-i)i for \z - i\ < —. 
i^o 

*The analytic continuation of the free resovent in the region \C\ < |, which avoids the introduction of the pole at 
= 0, can be constructed via the formula 



Jo 



The energy estimate \\dt /II2 < II/II2 together with the sharp Huyghens principle lead to the bound 

II ||2~>2 ^ e"^*''^, which is sufficient for analytic continuation. 
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where the coefficients are compact operators in L?. Indeed, in Lemma 17.71 (and its generahzation 
to n > 3 odd) we observed that Ffr{i — is analytic for |C| < |. Thus, the series in ()7.5U|) is 
absolutely convergent in the slit region Ui{TZi) for small z. Moreover, since the rate of decay for the 
free evolution e**'^ is known to be the argument from the proof of Proposition 17. 9| applied 

to the free evolution, yields that all Cj = for 1 < j < n — 2 odd. Finally, by the resolvent identity 
one has, see (|7.24|) 



Gs{z) = (l + Fe{z)E-'V] F,{z 



-1 



-1 



-1 



(7.51) = (l + F,{i)E-'Vj [l + {F,{z)-F,{i))E-'V(l + F,{i)E-'Vj J F,{z). 

This is justified, since by our assumption of iti being neither eigenvalue nor resonance of A the 
inverse (l+Fs{i)E-'^v) exists as a bounded operator on L^, cf. (|7.25|1 . Furthermore, the inverse 



in square brackets in 1)7. 51(1 exists as an convergent Neuman series. Inserting ()7.5U|) into (|7.51() 
reveals that all coefficients of the entire right-hand side in (|7.5H) corresponding to powers {z — i)^ 

n — 2 

with 1 < i < n-2 odd, have to vanish. Therefore, ((7^ holds with \\C{t)\\ < \T-i\~, and ifT^ 
thus holds with {t)~^ . 

The third place where the dimension becomes relevant is the argument allowing for removal 
in Theorem 17.41 As in the scalar case, a tedious calculation involving further Duhamel expansion 
would extend the 3-d argument to higher dimensions. 

8 Matrix charge transfer model 

Definition 8.1. By a matrix charge transfer model we mean a system 



.1) 



tp\t=o = V'o, 

where Vj are distinct vectors in M^, and Vj are matrix potentials of the form 

/ Ujix) -e'OAt,^)W,ix) 

where 6j{t, x) = {\vj p -|- a|)t -|- 2x • Vj + 7j, aj,jj G M, aj ^ 0. Furthermore, we require that each 
be admissible in the sense of Definition\7.'J^ and that it satisfy the stability condition ()7.11|) . 



In comparison to the simple definition of a charge transform model in the scalar case, Defini- 
tion |H3 might seem unnatural. However, it is natural by virtue of being the only Galilei invariant 
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definition. In our context we need to use the following vector- valued Galilei transform (which 
should explain why we are using (|7.5() rather than ((731)): 

(.2) 



where g^^y are as in Section 2. Clearly, the transformations Qv,y{'t) are isometrics on all spaces. 
Since in our case always y = 0, we set Q^{t) := Qvfi{t). By 1)2. 7|) . Q^it)^^ = G-^(t) in that case. 
There is now the following analogue of ()2.7() . ()2.8|) . In contrast to the scalar case, the transformation 
law of Lemma 18.21 also involves a modulation M{t). It should be compared with Definition 12.11 



Lemma 8.2. Let a E M and set 

A : = 



with real-valued U, W. Moreover, let v £ M^, 6{t, x) = (|w p + Q^)t + 2x • i; + 7, 7 € M, and define 

, . _ f -iA + [/(• - vt) -e*''(*'-^^*)VF(- - vt 
^ ' V e-^^(*'-^^*)VF(- - vt) iA - U{- - vt) 



Let S{t), 5(0) = Ld, denote the propagator of the system 

1 



dtS{t) + H{t)S{t) = 0. 



Finally, let 



(8.3) M{t) = Ma,,{t) = Q 
where uj(t) = o?t + 7. T/ien 

(8.4) 5(t) = g,^(t)-i7W(t)-ie-^*^7W(0)g,^(0). 
Proof. One has 

(8.5) iaiA^(t)g,^(t)5(t) = f "l'^ °. ^A^(t)g,^(t)5(t)+A^(i)ig^-(i)5(t)-A^(t)g„^(t)i/(t)5(t). 



Let p{t,x) = tl^ P + 2x • One now checks the following properties by differentiation: 

M{t)-g^{t) = ni^l' + ^-p \M{t)g^{t) 

I \ —-^Ivl +v-p J 

( --A + U -e*(''~^~'^)w \ 

(8.6) + _,|.|»^_ _)>,(,)e,(t). 

In view of our definitions, 6 — p — lu = 0. Since co = a^, the lemma follows by inserting 1)8. 6() 
into dSISI)- □ 
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We now return to the matrix charge transfer problem (jH-lj) . In order to state our main theorem, we 
need to impose an orthogonahty condition in the context of the charge transfer, analogous to the 
one used in the scalar case. To do so, let Pc{Hj) and Pb{Hj) be the projectors as in Lemma 17.31 
Abusing terminology somewhat, we refer to Ran(Pfe(iJj)) as the bound states of Hj. Since Pb 
and Pc are no longer orthogonal projections, we use "scattering states" instead of "asymptotically 
orthogonal to the bound states" in the following definition. 

Definition 8.3. Let U(t)ipo = ip{t, •) be the solution of (|8.1|) . We say that ipo is a scattering state 
relative to Hj if 

\\Pb{Hj,t)U{t)^po\\L2 ^0 ast^ +oo. 

Here 

(8.7) Pb{Hj,t) := {t)-^Mj{t)-^Pb{H,)Mj{t)g^^ (t) 
with Mj{t) = Maj,jj{t) as in (|8.3j) . 

The formula (|8.7|) is of course motivated by (|8.4j) . Clearly, Pb{Hj,t) is the projection onto the 
bound states of Hj that have been translated to the position of the matrix potential Vj{- — tvj). 
Equivalently, one can think of it as translating the solution of 1)8. 1|) from that position to the origin, 
projecting onto the bound states of Hj, and then translating back. 
We now formulate our decay estimate for matrix charge transfer models. 

Theorem 8.4. Consider the matrix charge transfer model as in Definition \8.1l Let U{t) denote 
the propagator of the equation (|8.1|l . Then for any initial data tpo £ n L"^ , which is a scattering 
state relative to each Hj in the sense of Definition \8.°A one has the decay estimates 

(8.8) ||f/(i)^o||L2+L- < (t)-i||V^o||LinL2. 

If in addition the matrix potentials Vj satisfy ||V^||li < oo for all j = then one also has the 

bound 

(8.9) ||f/(t)^o||L- < Itr^llV^ollLinL^. 

Remark 8.5. Theorem 18 . 41 also holds in the case of higher dimensions re > 3. In particular, we have 
the corresponding estimate 

||^(i)V'o||L°° < Itl'^Uoh^nL^- 

As in the case of scalar charge transfer models the proof of Theorem 18.41 mainlv relies on the decay 
estimates for the corresponding matrix problems Hj with a single time-independent potential Vj{x) 

(8.10) ||e^*^^P.(/7,)/llL^+L°o < {tr'^WfWL^nL- 

We have established such estimates for the admissible Hamiltonians Hj in Theorem 17.41 for re = 3. 
In section 7.2 we have also observed that our method for deriving 1)8. 10() can be generalized to treat 
the case of an arbitrary odd dimension n > 3. While our method leaves the even dimensional case 
open at the moment, one can, in principle, substitute Cuccagna's dispersive estimates |Cuj instead. 
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Remark 8.6. The conclusions of Theorem 18.41 are also valid for perturbed matrix charge transfer 
Hamiltonians. This refers to equations of the type 

Tp\t=0 = V'O, 

where the charge transfer part is as in Definition 18. II and the perturbation satisfies 

sup \\Vo{t, Ollinoo < £• 

i 

See the scalar case in Section [3.71 for an exact formulation. 

As in the scalar case, (|8.8|) is proved by means of a bootstrap argument with the same assump- 
tion Since the argument is basically identical with the scalar case, we do not write it out 
in full detail. As in the scalar case, the estimate ()8.9() is a consequence of the L2 + bound 
(|8.8|) . This follows by combining the machinery developed in Proposition 15.11 with the cancellation 
property for a matrix problem with a single potential established in the proof of the estimate 1)7. 13() 
in Theorem 17.41 

The fact that systems have generalized eigenspaces rather than just eigenspaces leads to some minor 
changes from the scalar case. But this really only affects the proof of Proposition 13.11 We show 
below that the statement still remains the same. 

In the following we shall assume that the number of potentials is z/ = 2 and that the velocities 
are iTi = 0,V2 = (1, 0, ... 0) = el. This can be done without loss of generality. 

8.1 Bound states 

We now estimate the rate of convergence of the projections onto the bound states of solutions which 
evolve from scattering states. 

Proposition 8.7. Let 'ip{t,x) = (U{t)i{jQ){x) be a solution of 1)8. 1(1 where ipQ is a scattering state 
relative to Hi and H2 in the sense of Definition \8.f\ Then 

||Pfe(i/l,t)C/(t)Vi'o||L2 + \\P^{H2,t)U{t)^4L^ < e-"*||V^o||L2 

for some a > 0. 

Proof. By symmetry it suffices to prove the bound on the first part (one can switch the roles of Vi 
and V2 by means of a Galilei transform and a modulation). In view of Lemma l7.3l one decomposes 

M 

(8.11) ${t) :=Ma,,jAt)Umo = J2fj{t)+Mt,-) 

j=0 

relative to Hi so that (j)i{t,-) lies in the continuous subspace of Hi, i.e., Pc{Hi)(pi = cpi and 
Pb{Hi)(pi = for all times. Furthermore, fj{t) G Lj where Lj are the spaces from Lemma 17.31 for 
the operator Hi. By assumption, 

M 

(8.12) ^ |/5(t)|2 ^ as t ^ 00. 

j=0 
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One checks that (p satisfies the equation 

(8.13) -dt$+Hi$+V2{--eit)^=0 

where V2 = Mar,-f^{t)V2Ma^,'f^{ty^ ■ Substituting l|H?TT|) into (|5Tn|l yields 

-dt$+Hi$+V2{--tei)$ 
I 

^ I ■ I ■ ^ 

(8.14) =^-f'j + -Mt) + +Y1 + ^i'?! + ^2(- - iel)< 

j=0 * ^ j=0 

Since Pc{Hi)(f)i = 0i and Ran(Pc(-ffi)) is closed by assumption, one has 

Hi$i = P^{Hi)Hi$i, dt$i = Pc{Hi)dt$i. 

In particular 

Pfe(/7i) + =0. 

Therefore, applying Pi,{Hi) to (|8.14p yields 



l^fjit) + ^i/5(0] = m := -PbiHi)V2i- - tel 



- 1 

j=0 



In view of the explicit expression (|7.8|) for P},{Hi) one has |5(t)| < e Moreover, decomposing 



^(t) = J2fLo9j{'^) with 5j(t) G ii, one also has |5j(t)| ^ e ^* for all < j < Af. Hence, 



(8.15) -fo + Hifo = go 

-Jj+^jfj = 9j for l<j<M 

with exponentially decaying right-hand sides. It follows from ()8.12|) that \fj{t)\ < e^^* for 1 < j < 
M. Moreover, applying Hi to (|8.15|) one obtains 

Since the right-hand side decays exponentially, and since -ffi/o(t) — > as t — > 00 by our assumption 
of starting from a scattering state, it follows that Hifo^t) decays exponentially, and therefore 
also fo{t). □ 

9 Generalized decay estimates for the charge transfer model 

Consider the time-dependent matrix charge transfer problem 

idtif+H{a,t)iJ=F 
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where the matrix charge transfer Hamiltonian H{a, t) is of the form 



^ 2 / 



3-- 

where Vj are distinct vectors in M^, and Vj are matrix potentials of the form 



Vj{t,x) 



e 



Uj{x) -e'^^^^'""^ Wj{x) 

-Uj{x) 



where 9j{t, x) = {\vj \ + aj)t + 2x ■ Vj + 7j, aj,^j G M, aj / 0. Our goal is to extend the dispersive 
estimate 

(9.1) WmWi^+L^ < (1 + i)"^ (llV^olliinL^ + |||F|||| + i?) 
with 

(9.2) |||F||| :=sup /V(^)llLidT + (l + t)t+i||F(t,.)||L2 

t>0 Jo 

to the corresponding estimates for the derivatives of il^{t). The estimate (|9.ip holds only for the 
solutions ip{t) which are scattering states, i.e., for ip obeying the a priori condition that 

\\PbiH,,t)i;{t)\\L2<B{i + ty^ 

for all J = 1, ..jV. Our first lemma shows that the functions 

4(t):=vV(t), keZ+ 

are scattering states as well. 

Lemma 9.1. The functions ipkit) obey the estimates 

(9.3) \\Pb{Hj,t)MmL^<Ck{i+tr'^ 



Proof. Let ff{x) be an arbitrary C°° exponentially localized function. Then for any y £ 

Mt,x)-ff{x-y)dx = {-l)'' [ V^(t,x)V'=7/(x-y)dx< ||V^(t)||i2+ioo||V^r/||iinL2 < (l + t)"^ 
Now recall that the projection 

Pb{Hj,t) := g,-yt)-'M,{tr^Pb{H,)Mjit)g,^yt) 

with the Pb{Hj is given exlicitly 

Pb{Hj)f = ^ CafSUaif, Vf3) 

where Ca/3 are given constants and Ua^vp are exponentially localized functions. The result now 
follows. □ 
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Proposition 9.2. The functions tj^k = V'^V satisfy the + L°° dispersive estimate 

k 

(9.4) IIvVWIIl^+z.^ < (1 + (llvVollLinL^ + II|V'F||| + b) 



i=0 

7k^ 



Proof. We have already shown that V is a scattering state. Moreover, differentiating the equation 
k times we obtain 



fc-i 

idtV'^iJ + Hit, cj) V V =Fk-=Y^ Gi{t, x) V V + V^F 

£=0 



where Ge{t, x) are smooth exponentially localized potentials uniformly bounded in time. Therefore 
V*^ is a scattering state solving an inhomogeneous charge transfer problem. Using the estimate 
(|9.H1 we then have 

(9.5) IIvVWIIl^+loo < (1 +t)-t(||vVo||LW + \m{r)\\\+B) 

We use that for any p G [1,2] 

||G,(t,x)vVl|Lp < ||vV||l2+l- 
Proceeding by induction on k we conclude that for any £ < k 

f \\G,{t,x)V'nL^ < Al + T)-t(ir^ (||V-V^o||LinL2 + ll|V"^F|||) 
•^0 JO ^^0^ ^ 

i 



m=0 

and that 



m=0 



The result now follows from 1)9. 5() and the inequality 



fc-i 



\Fk{T)\\\ < |||V'=F|||+^|||G,(t,x)V^ 



1=0 



□ 



Consider the following Banach spaces X and y of functions of {t,x): 

s 

(9.6) IIV^IU^ = sup(||V(t,-)||H»+(l+i)^Ell^'^(*'-)lli^+^' 

(9.7) ||Fb^ = sup^( r||V'i^(r,-)llL^dT + (l + t)i+^||V'^F(t,.)||L^ 



These spaces are relevant for our companion paper |RSSj . We can summarize our estimates for the 
charge transfer model in the following proposition. 
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Proposition 9.3. Let ip be a solution of the matrix charge transfer problem 



idttp + H{t, a)ip = F 



satisfying the condition that for every j = 1, 



(9.8) 



\\Pb{Hj{a,t))^\\L2<B{l + tr^ 



Then for any integer s >0 



s 



(9.9) 



k=0 
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